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I.  INTROUJCTIOh 


Tnis  is  a  final  report  on  a  three-year  fellowship  program  in  hypersonic 
ae -ody na~ i cs  sporsorec  at  the  University  of  Maryland  by  the  Army  Research 
Office.  Tnis  procram  was  initiated  in  '-'.arch  19o3,  and  was  an  insightful 
precurser  to  the  massive  activity  in  hypersonic  aerodynamics  now  sweeping  the 
country.  This  ARO  fellowship  program  was  characterized  by  the  following: 

(1)  It  was  the  first  fellowship  program  in  hypersonic 
aerodynamics  in  an,  university. 

(a,  It  helped  to  spa'v  c  researcn  activity  in  hypersonic 
ae-c:..  na  -  i  cs  at  tr-e  v-5  versity  o'  Ma -ylard--an  activity  that  now 
e~cc%  asses  1?  f  .ill -time  graojate  stjdents. 

(3;  It  sponsored  two  exceptional  Rr.O.  students--both  straight  A 
students--*’ tn  fellowships  in  hypersonic  aerodynamics  for  three 
years. 

[i\  Qre  0<-  the-ne  students,  Hr.  Kevin  Bowcutt,  graduated  in  August 
19.-M.  Me  has  since  gone  on  to  distinguish  himself  as  a  leader  in 
the  hat i on  a  1  Aerospace  Plane  procra~  at  Rockwell  International. 

(S'  Tne  othe-  student  is  Griffin  Corpening,  who  will  graduate  in 
May  19co.  M' .  Corpening  is  presently  working  at  the  Jonns  Hopkins 
App'ie:  R-ysics  Laot-at  :>rj  in  tnei  r  hypersonic  vehicle  program. 

I r  light  o*  tne  amove,  tne  ARC  Fellowship  Program  in  Hypersonic  Aerodynamics  has 
rac  a  synergestic  efiect,  not  only  at  Maryland,  but  on  the  national  hypersonic 
p'-ogra"’  in  the  United  States  far  beyond  the  initial  expectations  for  the 
nrogra~. 

In  the  following,  the  research  conducted  under  this  fellowship  program  is 
'w  an  zed. 


1 1  .  \!'.$C?"S  QPTIVIZEl  HYPERSONIC  U'AVER I DCR5 

Here  the  research  is  summarized  by  two  AIAA  papers,  reproduced  in  total  o 
the  followi no  page-.  One  paper,  by  Bowcutt  and  Anderson,  represents  work 
totally  sponsored  by  ARO,  and  is  derived  from  the  Ph.D.  dissertation  of  Kevin 
Bowcutt.  Tnis  paper  is  self-explanatory.  Tne  second  paper,  by  Corda  and 
Anderson,  represents  a  deri vation  of  the  ARO  sponsored  work,  and  is  included 
here  to  demonstrate  the  synergestic  effect  of  this  grant. 
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VISCOUS  OPTIMIZE  HYPERSONIC  WAVERIDERS 
DESIGNED  FROM  AXISYMMETRIC  FLOW  FIELDS 
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ABSTRACT 

A  series  of  "viscous  optimized"  waveriders  is 
designed  from  general  axi symmetric  flow  fields, 
otner  than  conical  --  most  notibly,  flow  fields 
over  minimum  drag  bodies.  Tnese  configurations 
represent  the  next  logical  step  in  an  ongoing 
program  of  research  on  the  design  of  hypersonic 
waveriders  at  the  University  of  Maryland.  In  pre¬ 
vious  work  at  Maryland,  inviscid  conical  flow  was 
used  exclusively  in  the  design  of  waveriders. 

These  conical  flow  waveriders  predict  high  values 
of  lift/drag,  higher  than  the  values  based  on 
experience  from  other  hypersonic  configurations. 
The  present  work  allows  the  design  of  hypersonic 
waveriders  from  any  axisymmetric  flow  field  — 
although  the  present  results  focus  on  the  inviscid 
flows  over  cones  and  3/4  and  1/2  power-law  bodies. 
A  space  marching,  finite-difference  code  is  used 
to  generate  the  axisymmetric  flow  field.  The 
undersurface  of  the  vehicle  is  carved  out  as  a 
stream  surface  of  this  axisymmetric  flow  field, 
whereas  the  upper  surface  is  assumed  to  be  a 
freestream  surface.  The  detailed  viscous  effects 
are  included  within  the  optimization  process, 
using  a  simple  reference  temperature  method. 
Results  obtained  using  the  reference  temperature 
method  are  within  101  of  results  obtained  using  a 
more  complex  integral  boundary  layer  method,  even 
at  high  hypersonic  Mach  numbers.  Boundary  layer 
transition  is  predicted  using  a  correlation  of  the 
local  transition  Reynolds  number  with  the  Mach 
nunber.  A  non-linear  simplex  method  is  used  to 
optimize  the  waveriders  for  either  maximum  L/D  or 
minimum  total  drag. 

INTRODUCTION 

In  recent  years,  research  involving  all 
aspects  of  hypersonic  flight  has  seen  a  rebirth  in 
the  United  States.  An  extended  survey  of  hyper¬ 
sonic  aerodynamic  research  is  given  in  Reference 

1.  During  the  long  hiatus  that  hypersonic 
research  lay  dormant,  powerful  new  engineering 
tools  have  become  available  --  i.e.  supercomputers 
and  computational  fluid  dynamics  (CFU).  In  fact, 
vehicle  concepts  such  as  the  National  Aerospace 
Plane  (NASP),  aero-assi sted  orbital  transfer 
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vehicles  (AOTV's),  the  hypersonic  transport  (the 
"Orient  Express"),  and  hypersonic  missiles,  place 
much  of  the  design  burden  upon  CFO.  In  addition, 
advanced  computers  and  computational  techniques 
allow  us  to  take  a  “new"  look  at  some  "old"  hyper¬ 
sonic  aerodynamic  concepts.  The  present  paper  is 
part  of  a  continuing  hypersonics  research  program, 
at  the  University  of  Maryland,  that  does  just 
that.  Using  modern  computational  techniques,  the 
"old"  concept  of  waveriders  is  revisited  to  pro¬ 
duce  a  "new"  class  of  advanced  hypersonic  lifting 
configurations.  This  CFD  approach  to  the 
waverider  concept  has  lead  to  a  design  tool  with 
the  following  "new"  advantages: 

1.  The  waverider  configurations  can  be  numeri¬ 
cally  optimized  for  almost  any  figure  of 
merit  --  maximum  L/D  and  minimum  total 
drag,  in  the  present  study. 

2.  The  detailed  viscous  effects  are  included 
within  the  optimization  process. 

3.  The  waverider  configurations  can  be  carved 
out  as  stream  surfaces  of  almost  any  axi¬ 
symmetric  flow  field  --  the  inviscid  flow 
over  cones  and  3/4  and  1/2  power-law  bodies 
are  used  in  the  present  investigation. 

4.  Even  with  the  above  detailed  effects  and 
versatility,  the  waverider  optimization 
procedure  is  computational  efficient, 
allowing  parametric  design  studies  to  be 
made. 

The  advantages  listed  above  will  be  elaborated 
upon  in  the  paragraphs  to  come. 

This  paper  is  a  sequel  to  Reference  2.  In 
References  2  and  3,  Bowcutt  and  Anderson 
established  a  new  class  of  hypersonic  vehicles 
that  broke  the  "L/D  barrier"  as  depicted  in 
Figure  1.  This  barrier  represents  a  general 
empirical  correlation  for  (l/D)max  ,  as  the 
freestream  Mach  number,  Increases  across  the 
supersonic  and  hypersonic  regimes,  based  on  actual 
flight  vehicle  experience.  It  is  given  by  the 
relation 

4 ( M  +  3) 


(obtained  from  Reference  4),  and  is  shown  by  the 
solid  line  in  Figure  1.  A  number  of  data  points 


for  optimized  waveriders  generated  by  Bowcutt  and 
Corda  are  shown  along  with  various  other  hyper¬ 
sonic  conf igurations  at  various  Reynolds  numbers 
and  Mach  numbers.  (The  points  on  this  graph 
correspond  to  different  references  as  given  in 
References  5  and  6.)  The  waveriders  generated  by 
8owcutt  are  seen  to  follow  the  relation  given  by 

S(M  +  2) 

=  - pi - 

00 

shown  by  the  dashed  curve  in  Figure  1.  The  pre¬ 
sent  waverider  results,  as  represented  by  the 
solid  circles  in  Figure  1,  will  be  discussed  in 
subsequent  sections. 

To  fully  appreciate  the  logical  sequence  of  the 
present  waverider  research  and  how  it  represents  a 
contribution  to  the  state-of-the-art ,  let  us 
review  the  fundamental  concepts  of  waveriders.  In 
1959,  Nonweiler7  proposed  the  idea  of  constructing 
a  three-dimensional  hypersonic  vehicle  from  a 
known  flow  field.  Nonweiler  chose  the  flow  field 
behind  a  planar  oblique  shock  wave  to  generate  a 
class  of  vehicles  with  a  caret-shaped  cross- 
section  and  a  delta  planform.  The  construction  of 
a  caret  wing,  as  they  are  called,  is  shown  in 
Figure  2  ( from  Reference  8).  When  flying  at  its 
design  Mach  number,  the  body  appears  to  be  riding 
on  top  of  the  attached  shock  wave,  and  hence  is 
dubbed  a  "waverider11.  There  is  no  flow  spillage 
from  the  lower  to  the  upper  surface  of  the  vehicle 
when  it  is  flying  at  its  design  Mach  number 
because  the  shock  wave  is  attached  to  the  leading 
edges.  This  containment  of  the  flow  beneath  the 
vehicle  results  in  high  pressure  being  exerted  on 
the  lower  surface  --  the  high  pressure  behind  the 
two-dimensional  planar  shock  wave.  This  high 
pressure  leads  to  high  lift.  It  can  be  stated 
that,  in  general,  when  compared  at  the  same  lift 
coefficient,  caret  wings  have  higher  values  of  I/O 
than  other  winged  hypersonic  vehicles.  The  aero¬ 
dynamic  benefits  of  waveriders  are  listed  in 
Reference  1,  and  are  discussed  in  detail  in 
References  5,  6,  and  9. 

The  idea  of  generating  a  three-dimensional 
hypersonic  vehicle  from  a  known  flow  field  can  be 
extended  from  Nonweiler's  construction  (which  uses 
a  simple  two-dimensional  wedge)  to  more  complex 
flows.  For  example,  a  conical  flow  can  be  used, 
as  shown  in  Figure  3  (from  Reference  8).  The 
resulting  class  of  "conical-flow"  waveriders  was 
first  investigated  by  the  British  (see  Reference 
10).  Rasmussen  et  al.  »  ,  Cole  and  Zien  , 

and  Kim  et  al.  **  represent  further  variations  — 
designing  waveriders  from  flows  over  cones  and 
elliptic  cylinders  using  hypersonic  small  distur- 


procedure,  a  configuration  can  be  found  which 
balances  the  desire  for  aerodynamic  efficiency 
(high  L/D)  with  the  large  wetted  area  (hence  high 
friction  drag)  attendent  with  this  high  L/D. 

The  present  work  takes  the  next  logical  step, 
going  beyond  conical  flow  fields,  to  generate 
"viscous  optimized"  waveriders  from  general 
axisymmetric  flow  fields  —  most  notibly,  flow 
over  minimum  drag  bodies.  Again,  the  detailed 
viscous  effects  are  included  within  the  optimiza¬ 
tion  process.  The  rationale  behind  this  direction 
is  as  follows.  One  might  intuitively  suspect  that 
the  flow  over  a  minimum  drag  power-law  body  might 
lead  to  waveriders  of  lower  drag  and  higher  L/D 
than  those  derived  from  conical  flow,  by  nature  of 
the  generating  flow  field  itself.  The  need  for 
minimum  drag  and  high  L/D  are  of  paramount  impor¬ 
tance  for  vehicles  (such  as  the  aerospace  plane) 
that  must  be  able  to  accelerate  and  cruise  as 
efficiently  as  possible.  One  might  also  expect  a 
better  volume  distribution  (meaning  more  volume  in 
the  nose  of  the  vehicle  for  packaging  of  payload 
or  avionics)  with  a  power-law  derived  waverider 
than  with  a  conically  derived  waverider.  In  par¬ 
ticular,  the  present  results  show  waveriders 
designed  from  the  flow  over  cones  and  1/2  and  3/4 
power-law  bodies.  Also,  the  configurations  are 
numerically  optimized,  using  the  npn-linear 
simplex  method  of  Nelder  and  Mead1  ,  for  maximum 
L/D  and  for  minimum  total  drag. 

ANALYSIS 

The  analysis  and  optimization  of  waveriders 
from  general  axisymmetric  flow  fields  require 
distinct  capabilities.  First,  one  must  be  able  to 
generate  a  "known"  axisymmetric  flow  field.  Then 
one  must  analyze  the  waverider  generated  from  the 
flow  field  which  involves:  1)  the  evaluation  of 
the  aerodynamic  coefficients  (e.g.  Cl,  Cg,  and 
Cm)  and  2)  the  determination  of  the  parameters 
which  describe  the  waverider  geometry.  Finally, 
waverider  optimization  requires  finding  the  set  of 
parameters  that  yield  the  desired  figure  of  merit 
(e.g.  maximum  L/D  or  minimum  Cp  ).  The  breakdown 
of  the  waverider  design  process  is  shown  in  Figure 
4.  The  contents  of  each  box  will  be  discussed  in 
the  following  paragraphs.  A  more  detailed 
description  of  the  design  process  is  given  in 
Reference  6. 

A.  Generation  of  Inviscid  Flow  Field 

The  calculation  of  the  "known"  axisymmetric 
flow  field  will  be  discussed  first.  Separate 
techniques  are  used  for  the  case  of  1)  flow  over  a 
right  circular  cone  and  2)  flow  over  a  general 
axisymmetric  body,  as  described  below. 


bance  theory.  Further,  these  investigators  used 
the  calculus  of  variations  to  search  fo£  ^ptimum 
waverider  shapes.  Bowcutt  and  Anderson  > 
generated  a  series  of  conical-flow  waveriders  that 
were  numerically  optimized  for  maximum  L/0. 

Unlike  previous  Investigators,  they  are  the  first 
to  include  the  detailed  viscous  effects  (including 
boundary  layer  transition)  within  the  optimization 
process.  This  helps  to  alleviate  a  fundamental 
drawback  of  waveriders  --  the  large  friction  drag 
associated  with  a  waverider's  characteristical ly 
large  wetted  area.  With  skin  friction  effects 
fully  accounted  for  within  the  optimization 


A.l  Conical  Flow 

For  the  case  of  conical  flow,  the  ordinary 
differential  equation  know  as  the  Taylor-Maccol 1 
equation  is  solved  using  a  standard  fourth-order 
Runge-Kutta  method  as  given  in  Reference  16  and 
shown  below. 


RiJI 


1 


?*v3 


a 


ft 

I 


>*V<rC 


•1‘Vr  *  W^^r  +  ^cotB+  -T1 

1  dh 

dV  dV  dV  dJV 

‘  if  -Vr  if  *  IT  ~f  1=0  (1) 

de 

In  equation  (1),  Vr  is  the  nondimensi onal  com¬ 
ponent  of  flow  velocity  along  a  conical  ray,  «  is 
the  angle  of  the  ray  referenced  to  the  cone  axis, 
and  y  is  the  ratio  of  specific  heats. 

A2.  General  Asymmetric  Flow  Field 

For  an  irtviscid  ax i symmetric  flow  field  in 
general,  the  governing  partial  differential 
equations  are  the  Euler  equations  given  helow 


G  -  p  +  pw 


H  *  y  cvw 


whe^e  p  is  the  density,  w  is  the  axial  ( z }  com¬ 
ponent  of  the  flow  velocity,  v  is  the  transverse 
(r)  component  of  the  flow  velocity,  and  p  is  the 
pressure. 

These  equations  are  solved  numerically  using 
MacCormack's  explicit,  space-marcing,  finite- 
difference  scheme.  In  this  technique,  the  flow 
field  of  interest  is  discretized  into  a  network  of 
grid  points  in  the  axial  (z)  and  radial  (r)  direc¬ 
tions.  Starting  from  an  initial  data  plane,  where 
the  flow  properties  are  assumed  known,  the  flow 
field  solution  is  calculated  by  marching 
downstream  in  steps  of  Az.  The  shock-fitted 
approach  is  used  where  the  solution  is  bounded  by 
the  surface  of  the  body  and  the  shock  wave  pro¬ 
duced  by  the  body.  Applied  to  the  Euler 
equations,  the  scheme  can  be  written  as 
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where  A z  is  the  marching  increment  in  the  axial 
(z)  direction,  Ar  is  the  distance  between  grid 
points  in  the  radial  (r)  direction,  n  is  the  grid 
point  index  in  the  axial  direction,  and  j  is  the 
grid  point  index  in  the  radial  direction.  To 
obtain  a  stable  solution,  the  downstream  marching 
is  limited  by  the  CFL  criterion 


az  =  C 


Ar 

tan(e  +  u) 


where  C  is  a  constant  between  0  and  1.  A  value  of 
C  equal  to  0.95  is  used  in  all  of  the  calculations 
presented  in  this  paper. 

As  mentioned  earlier,  the  flow  field  solution 
is  started  from  a  known  data  line.  For  this 
reason,  the  body  of  interest  is  assumed  to  have  a 
conical  nosetip.  The  flow  field  behind  the  conical 
nosetip  is  solved  using  the  Taylor-Maccol 1  cone 
flow  solution  outlined  in  the  previous  section. 
This  nosetip  flow  field  is  then  matched  with  the 
body  of  interest  to  serve  as  an  initial  data  line 
for  the  space-marching  solution.  For  a  detailed 
account  of  HacCormack 1 s  space-marching  technique, 
including  the  application  of  the  boundary  con¬ 
ditions,  see  References  6  or  17. 
body  of  interest  to  serve  as  an  initial  data  line 
for  the  space-marching  solution.  For  a  detailed 
account  of  MacCormack ' s  space-marching  technique, 
inrluding  the  application  of  the  boundary  con¬ 
ditions,  see  References  6  or  17. 

B.  Generation  of  Leading  Edge  Shapes 

Once  the  flow  field  is  calculated,  waverider 
shapes  can  be  "carved"  out  of  the  flow  field  as 
shown  in  Figure  5.  As  is  seen  in  this  figure,  the 
leading  edge  curve  is  traced  back  along  the 
streamlines  of  the  flow  field  to  generate  the 
lower  surface  of  a  waverider.  Note  that  the  pro¬ 
jection  of  the  leading  edge  curve  in  the  cross- 
stream  (x-y)  plane  uniquely  defines  a  waverider 
geometry  for  a  given  flow  field.  The  leading  edge 
curves  are  generated  by  the  optimization  process 
as  will  be  discussed  shortly. 

C.  Streaml ine  Tracing 

As  noted  in  the  previous  section,  stream¬ 
lines  are  traced  from  the  leading  edge,  through 
the  "known"  flow  field,  to  create  the  waverider 
lower  surface.  The  manner  of  this  streamline 
tracing  will  now  be  described. 

From  the  space-marching  solution,  flow 
properties  are  known  at  distinct  axial  (z)  sta¬ 
tions.  Tracing  a  streamline  amounts  to  finding 
where  the  value  of  the  streamfunction,  ip ,  along 
each  z  data  line  is  equal  to  the  value  of  the 

streamfunciton  of  the  specified  streamline.  The 
streamfunction  along  each  data  line  can  be  calcu¬ 
lated  by  using  the  definition  of  the  streamfunc¬ 
tion  for  an  axially  symmetric  flow. 


m 


Vi  =  *j  + [b] 

where  •')  is  the  value  of  tne  streamfunction,  j  is  a 
radial  (r)  index  for  data  points  along  a  z 
station,  and  n  is  the  mass  flow  in  the  annulus 
bounded  by  the  points  j  and  j+1.  8y  setting  the 
value  of  the  streamfunction  at  the  body  surface 
equal  to  zero,  the  value  of  the  streamfunction  can 
be  calculated  along  each  z  data  line. 

0.  Freestream  tipper  Surface 

The  generation  of  the  waverider  lower  sur¬ 
face  by  streamline  tracing  has  been  described  in 
the  previous  sections.  The  upper  surface  is 
created  by  simply  following  the  freestream  back 
tnrough  the  given  leading  edge  curve  to  the  base 
of  the  waverider.  The  pressure  on  this  freestream 
upper  surface  is  the  freestream  pressure,  p^. 

E.  Skin  Friction  Calculations 

The  skin  friction  distribution  along  the 
streamlines  that  form  the  waverider  is  calcinated 
us’og  the  reference  temperature  method  of 
Eckert'9.  In  the  reference  temperature  method, 
approximate  formulas  are  used  to  predict  the  skin 
friction,  with  the  physical  properties  evaluated 
at  an  appropriate  reference  temperature.  For  a 
flat  plate  in  laminar  flow,  the  local  skin  fric¬ 
tion  coefficient  is  given  by 

cf  =  0.664  (*ex)-l/2  (T7T.)(u*1)/2  (6) 

Rex  is  the  local  Reynolds  number  defined  as 
p  V  x 


where  o«  is  the  freestream  density,  V„  is  the 
freestream  velocity,  x  is  the  local  distance  from 
the  leading  edge  of  the  plate,  and  u„  is  the 
freestream  value  of  the  viscosity. 

Also  in  equation  (6),  T*  is  the  reference  tem¬ 
perature  defined  as 

(T 7T J  =  1  +  0.032  M2  *  0.58  (T^/T,  -  1)  (8) 


where  H»  is  the  freestream  Mach  number  and  Tw  is 
the  wall  temperature.  Finally,  w  in  equation  (6) 
is  the  exponent  of  an  assumed  exponential 
variation  of  u ,  namely, 

(u’/uj  -  (T7TJ“  (?: 

A  value  of  u  =  0.75  is  used  in  the  present  study. 

The  flat  plate  skin  friction  coefficient  for 
turbulent  flow  is  given  by 

.  0.0592  (10a) 


p'  V  x 


In  equation  (10b),  .  1  and  are  evaluated 
at  the  reference  temperature  T7 

This  type  of  skin  friction  analysis  is  much 
simpler,  in  concept  and  use,  than  integral  boun¬ 
dary  layer  methods  (such  as  used  in  Reference  2). 
Results  obtained  using  the  reference  temperature 
method  were  within  10%  of  results  using  a  more 
complex  integral  boundary  layer  method  (Reference 
2),  even  at  high  hypersonic  Mach  numbers.  Also, 
the  computation  time  required  by  the  reference 
temperature  method  is  very  small  when  compared 
with  an  integral  boundary  layer  method. 

Boundary  layer  transition  is  predicted  using  a 
correlation  of  the  local  transition  Reynolds 
number,  Rex  t  with  the  local  edge  Mach  number, 

Me,  as  follows 

log.-Re  =  6.421  e^1'209*10  “Me  J  (11) 
•’10  xt 

This  correlation  is  the  same  as  used  in  the  pre¬ 
vious  work  by  Sowcutt  and  Anderson2’3  and  is  based 
on  the  experimental  data  of  OiCristina  9  for  tran¬ 
sition  on  sharp  cones  at  zero  degrees  angle  of 
attack.  This  correlation  is  used  due  to  the  lack 
of  better  methods  of  transition  prediction  in 
hypersonic  flows. 

F.  Aerodynamic  Forces  and  Moments 

The  aerodynamic  force  and  moment  coefficients 
are  calculated  by  numerically  integrating  the 
pressure  and  shear  stress  over  the  surface  of  the 
waverider.  Base  drag  is  not  calculated.  Details 
of  the  integration  for  the  forces  and  moments  can 
be  found  in  Reference  6. 

G,  Waverider  Optimization 

The  construction  and  aerodynamic  analysis  of  a 
single  waverider  configuration  is  shown  by  blocks 
A  through  F,  in  Figure  4,  and  just  described  in 
the  previous  paragraphs.  Now,  the  optimization  is 
performed  by  perturbing  the  shape  of  the  leading 
edge  curve  (which  corresponds  to  a  unique 
waverider  geometry),  until  a  configuration  is 
found  with  the  optimum  value  of  the  specified 
figure  of  merit  —  either  maximum  L/D  or  minimum 
total  drag,  in  the  present  analysis.  Th^non- 
1  inear  simplex  method  of  Nelder  and  Mead  is  used 
for  the  numerical  optimization.  Because  this  is  a 
zero  order  method,  only  function  evaluations  are 
needed  to  find  an  optimum  (no  derivative  infor¬ 
mation  is  needed).  These  function  evaluations  are 
the  determinations  of  the  figure  of  merit  (e.g. 

L/D  or  Cg)  for  each  configuration.  The  zero  order 
methods  are  the  easiest  to  implement  and  it  was 
shown  by  VanWie20  that  the  non-linear  simplex 
method  worked  as  well  as  other  higher-order 
methods  when  applied  to  a  similar  type  of  optimi¬ 
zation  problem  as  is  considered  here. 

The  non-linear  simplex  method  minimizes  a 
function  of  n  variables  by  comparing  values  of  the 
function  at  (n+1)  vertices  of  a  ’’simplex", 
replacing  the  vertex  with  the  highest  function 
value  with  another  point  determined  by  the  logic 
of  the  scheme,  and  then  "moving"  this  simplex, 
over  the  function  surface,  in  the  direction  of  the 
function  minimum.  Three  operations  are  used  by 


levels  are  used  to  reach  an  "optimium"  waverider 
configuration. 


the  simplex  in  its  search  for  the  function  minimun 
--  reflection,  expansion,  and  contraction.  These 
operations  are  graphically  illustrated  by  con¬ 
sidering  the  minimization  problem  shown  in  Figure 
6.  Here,  the  minimun  of  a  function  of  two 
variables  (X^  and  Xg)  is  sought,  where  the  contour 
map  defines  the  function.  A  simplex  composed  of 
n+l=3  vertices  (a  triangle)  moves  over  the  func¬ 
tion  surface  by  reflecting  or  "flip-flopping"  down 
the  function  valley,  expanding  if  possible  to 
speed  up  the  process,  and  finally  contracting 
around  the  function  minimum.  Constraints  in  the 
optimization  process  take  the  form  of  "barriers" 
or  boundaries  that  the  simplex  is  not  allowed  to 
cross.  The  nature  of  these  constraints  will  be 
discussed  in  a  later  section. 

In  the  present  worx  ,  the  function  to  be  mini¬ 
mized  is  the  negative  of  the  I'ft/drag  ratio  or 
the  value  of  the  total  drag  coefficient.  Fac"> 
function  evaluation  corresponds  to  a  unique 
waverider  configuration,  which  can  be  parame¬ 
terized  by  its  leading  edge  shape.  Therefore,  a 
leading  elge  curve  represents  a  single  vertex 
point  of  the  simplex.  The  leading  edge  is  spe¬ 
cified  by  five  points  in  the  crjss-stream  (x-y) 
plane,  where  a  cjbic  spline  is  used  to  define  a 
smooth,  continuous  curve.  It  should  be  noted  tnat 
calculations  are  performed  for  only  half  of  tne 
vehicle  because  of  its  longitudinal  symmetry.  One 
of  the  five  points  is  constrained  to  lie  on  the 
symmetry  plane  (y-axis)  of  the  waverider,  thus 
reducing  the  number  of  x,  y  coordinates,  needed  to 
define  the  leading  edge,  from  10  to  9.  The  number 
of  coordinate  variables  is  further  reduced  to  8  by 
the  constraint  that  the  last  leading  edge  point 
must  lie  on  the  shock  wave.  It  is  these  eight 
varialbes  that  are  perturbed  by  the  log’c  of  the 
optimization  routine  to  find  an  optimum  .  >'°rider. 

Since  the  optimization  problem  is  reduced  to 
finding  the  minimum  o*  a  function  defined  by  eight 
variables  ( n = 8 ) ,  the  simplex  for  this  problem  has 
n+1  or  nine  vertices  with  each  vertex  corre¬ 
sponding  to  a  set  of  coordinates  for  a  leading 
edge  curve.  To  initiate  the  simplex  method,  nine 
leading  edge  shapes  must  be  chosen  --  so-called 
basi s  leading  edge  shapes.  In  the  present  study, 
six  of  the  basis  leading  edge  shapes  are  poly¬ 
nomials  of  the  form 

yle  1  C1  *  C2xle  +  C3xle  +  C4xle  <12> 
and  three  of  the  basis  shapes  are  of  the  forms 

TT  X  . 

y1e  =  c5  +  c6  (1  -  cos  ~r~^  )  (13a) 

rX! 

yle  =  C7  +  c8  5in  *  - )  (13b) 

where  x ] e  and  y^e  are  the  x  and  y  coordinates  of 
the  leading  edge  and  rs  is  the  radius  of  the  shock 
wave  at  the  base  of  the  waverider.  A  set  of  basis 
leading  edge  shapes  is  shown  in  Figure  7.  The 
final  optimized  leading  edge  curve  is  also  shown 
as  the  bold  line  in  this  figure.  A  typical  opti¬ 
mization  history  is  shown  in  Figure  8.  In  all  of 
the  results  presented,  100  optimization  steps  or 


III.  RESULTS  AND  DISCUSSION 

The  presented  results  are  representative  of 
the  work  done  in  Reference  6.  The  following  areas 
will  be  covered:  (1)  validation  of  the  computer 
code  used  to  construct  and  aerodynamical ly  analyze 
a  waverider  (no  optimization);  (2)  a  presentation 
of  optimum  waveriders  designed  from  conical  and 
power-law  body  flow  fields,  over  a  wide  range  of 
hypersonic  Mach  numbers;  (3)  the  comparison  of 
waveriders  optimized  for  maximum  t/D  with  those 
optimized  for  minimun  total  drag;  (4)  the  effect 
of  using  an  all  laminar,  all  turbulent,  or  tran¬ 
siting  boundary  layer  in  the  waverider  optimiza¬ 
tion;  and  (5)  a  word  will  be  said  concerning  the 
computing  time  required. 

A.  Code  Val idation 

Because  of  the  unique  design  method  used  in 
the  present  analysis,  it  is  impossible  to  find 
experimental  data  to  confirm  the  computational 
results  presented.  However,  certain  point  com¬ 
parisons  can  be  made.  One  such  case  is  the  hyper¬ 
sonic  flow  over  a  half-cone  underneath  a  delta 
wing  configuration  at  zero  degrees  angle  of 
attack.  Experiments  by  Fetterman21  on  a  4  degree 
semi-apex  angle  half-cone  under  a  delta  wing  of  31 
degrees  sweepback  closely  approximate  tests  of  a 
waverider  configuration  at  a  Mach  number  of  6.86. 
Table  1  summarizes  the  results  of  this  comparison. 
The  wall  temperature  of  the  model  in  the  experi¬ 
ments  is  unknown,  so  calculations  are  made  for  a 
wall  temperature  equal  to  the  freestream  tem¬ 
perature  and  equal  to  the  total  temperature  of  the 
flow.  Also,  the  model  used  in  the  experiment  has 
a  1  degree  compression  on  the  upper  surface  of  the 
delta  wing  from  the  leading  edge  to  approximately 
the  midchord  position.  It  was  not  feasible  to 
model  this  compression  surface  so  calculations 
were  made  for  a  freestream  upper  surface  and  for  a 
1  degree  compression  on  the  entire  upper  surface. 
Table  1  shows  that  calculations  of  the  lift  coef¬ 
ficient,  C[_ ,  are  in  good  agreement  with  experi¬ 
ment.  Tne  calculated  values  of  the  drag 
coefficient,  Cq,  for  the  freestream  and  1  degree 
compression  upper  surface,  bracket  the  experimen¬ 
tal  value  as  would  be  expected.  This  same  trend 
follows  through  to  the  comparison  of  the  lift/drag 
ratio. 

Another  source  for  comparison  of  the  present 
results  is  the  previous  work  by  Bowcutt2,  for 
conical  flows  only.  This  is  possible  because  all 
of  the  computational  tools  used  in  the  present 
investigation  were  developed  independently  of  this 
previous  work.  Indeed,  for  the  case  of  conical 
flow  waveriders,  the  present  code  and  the  previous 
code  developed  by  Bowcutt  differ  by  a  few  percent 
for  most  parameters  of  Interest.  For  a  Mach 
number  equal  to  6  and  a  cone  angle  of  12  degrees, 
the  present  code  predicts  an  L/D  of  7.34  and 
Bowcutt's  code  predicts  an  L/D  of  7.73  —  a  dif¬ 
ference  of  about  5  percent.  For  a  higher  Mach 
number  of  14,  the  values  of  the  lift/drag  ratio 
are  6.06  and  5.87  from  Bowcutt's  code  and  the  pre¬ 
sent  code  respectively  --  a  difference  of  about  3 
percent.  This  comparison  also  helps  to  validate 
the  reference  temperature  method,  used  for  the 


viscous  flow  analysis  in  the  present  work,  since 
Rowcutt's  code  used  a  much  more  complex  integral 
boundary  layer  method.  The  comparison  is  espe¬ 
cially  encouraging  at  the  higher  Mach  numbers 
where  the  accuracy  of  the  reference  temperature 
method  may  be  suspect. 

3.  Waveriders  Designed  From  Cones  and 

Power-Law  Bodies 

Waveriders  are  optimized  at  specific  design 
points  (e.g.  Mach  number  and  altitude)  along  the 
possible  flight  trajectory  of  an  aerospace  plane. 
These  design  points  are  flight  Mach  numbers  of  4, 
6,  10,  14,  and  20  with  correspondi ng  altitudes  of 
80,000,  110,000,  125,000,  140,000,  and  175,000 
feet.  All  calculations  are  for  a  calorically 
perfect  gas. 

The  design  of  the  waveriders  are  optimized 
under  a  set  of  geometric  constraints  that  are 
easily  specified  by  the  program  user.  These 
constraints  include  the  following: 

1.  A  constraint  on  the  lower  limit  of  the 
slenderness  ratio  (reciprocal  of  the  fine¬ 
ness  ratio)  of  the  waverider--equal  to 
0.075  times  the  length  of  the  vehicle,  in 
the  present  results. 

2.  Constraints  on  the  box  size  (semi-span-to- 
length  ratio)  of  the  waveri der--the  lower 
limit  is  0.1  and  the  upper  limit  is  0.4,  in 
the  present  study. 

3.  A  constraint  on  the  minimum  acceptable 
total  volume  of  the  vehicle  --  not  an 
active  constraint  for  the  present  results. 

A  study  was  made  to  assess  the  effect  of  these 
constraints  on  the  design  process  but  will  not  be 
discussed  here.  Reference  6  can  be  consulted  for 
complete  detai Is. 

The  present  results  are  for  waveriders 
designed  from  flow  over  cones  and  1/2  and  3/4 
power-law  bodies.  Typical  body  shapes  are 
illustrated  in  Figure  9  for  a  Mach  number  of  10 
and  a  base  radius-to-length  ratio,  rbase/'- , 

0.15.  For  a  given  Mach  number  and  body  type,  the 
design  strategy  is  as  follows.  A  fineness  ratio 
of  the  flow  generating  body  is  selected  and  the 
flow  field  corresponding  to  this  body  is 
calculated  using  either  the  Tayl or-Maccol 1 
solution  (if  the  body  is  a  cone)  or  the 
space-marching  technique  (if  the  body  is  a  power- 
law  body).  The  optimum  waverider  is  then  found 
for  this  flow  field.  Additional  fineness  ratios 
are  selected  and  the  optimum  waverider 
correspondi ng  to  each  of  these  flow  fields  is 
obtained.  Figure  10  shows  a  series  of  waveriders 
optimized  for  maximum  I/O  for  a  Mach  number  of  6. 
Each  one  of  the  data  points  corresponds  to  a 
waverider  that  is  optimized  for  a  given  fineness 
ratio  of  the  flow  field  generating  body,  although 
the  lift/drag  ratio  is  plotted  versus  the  fineness 
ratio  of  the  waverider  itself.  From  Figure  10,  it 
is  seen  that  there  is  an  “optimum  of  the  optimums" 
for  each  of  the  body  types.  In  general,  it  is 
also  evident  that,  for  a  given  waverider  fineness 
ratio,  the  waveriders  designed  from  the  1/2 
power-law  body  have  the  highest  L/0,  followed  by 


those  designed  from  the  3/4  power-law  body  and 
those  designed  from  the  cone. 

All  of  the  remaining  results  in  this  section 
are  for  waveriders  that  are  optimized  for  maximum 
L/0.  Figure  11  shows  a  plot  of  the  lift/drag 
ratio  versus  the  Mach  number  for  waveriders  that 
are  the  “optimum  of  the  optimums".  In  general, 
the  power-law  derived  waveriders  have  slightly 
higher  values  of  L/0  than  the  conically  derived 
waveriders  over  the  Mach  number  range  shown.  The 
waveriders  designed  from  1/2  power-law  bodies  are 
seen  to  be  the  best  overall,  in  terms  of  high  L/D. 
Also  note  that,  for  all  the  body  types,  the  value 
of  L/0  decreases  then  increases  with  increasing 
Mach  number.  The  waverider  configurations 
designed  from  the  1/2  power-law  body  and  the  cone, 
for  Mach  numbers  of  4  and  20,  are  shown  in  Figures 
12  and  13. 

Figure  14  shows  a  plot  of  the  lift/drag  ratio 
versus  the  lift  coefficient,  C|_,  for  the  optimum 
waveriders  discussed  in  the  previous  paragraph. 
Note  that  decreasing  C|_  in  this  figure  corresponds 
to  an  increasing  Mach  number.  Here,  the  power-law 
body  derived  waveriders  generally  predict  higher 
values  of  L/D,  for  a  given  value  of  C|_ ,  than  the 
cone-derived  waveriders.  Again,  note  the  decrease 
then  increase  in  L/D  as  the  value  of  C(_  decreases 
(Mach  number  increasing). 

Figure  15  shows  L/D  plotted  against  a  Mach 
number  times  the  fineness  ratio  of  the  flow 
generating  body.  This  figure  shows  a  reversal  in 
trends,  in  terms  of  which  flow  field  yields  the 
waverider  with  the  highest  L/D.  The  power-law 
bodies  are  better  at  the  lower  values  of  M,(b/lb), 
and  the  cones  are  better  at  higher  values.  Also, 
it  is  seen  that  the  shapes  of  the  curves  in  Figure 
15  are  similar  for  each  type  of  flow  field.  This 
tends  to  suggest  that  the  quantity  M»(b/lb)  is  a 
valid  hypersonic  similarity  parameter  for  the  type 
of  waverider  design  process  presented  in  this 
paper. 

Returning  to  figure  1,  note  the  symbols 
denoted  by  solid  circles.  These  symbols 
correspond  to  the  waveriders  with  the  highest  L/0, 
designed  from  the  flow  over  1/2  power-law  bodies, 
at  the  various  Mach  numbers.  Note  that  they  also 
break  the  "L/0  barrier"  as  described  earlier.  As 
described  earlier,  the  data  for  L/D  displays  a 
“bucket"  type  behavior  with  increasing  Mach 
number.  The  decrease  then  increase  in  the  value 
of  L/D,  as  Mach  number  increases,  is  a  Reynolds 
number  effect  as  discussed  in  Reference  2, 

The  total  drag  coefficient,  Cg,  is  plotted 
against  Mach  number  for  the  "optimum  of  the 
optimums"  waveriders  in  Figure  16.  The  value  of 
drag  coefficient  is  seen  to  decrease  with 
increasing  Mach  number  for  all  of  the  flow  types. 
Note  that,  except  for  the  Mach  6  case,  the 
cone-flow  derived  waveriders  have  the  slightly 
higher  values  of  drag  coefficient.  It  is 
interesting  to  note,  in  Figure  17,  the  breakdown 
of  drag,  into  wave  and  skin  friction  drag,  for  the 
3/4  power-law  body  flow.  It  is  seen  that  the  skin 
friction  drag  and  wave  drag  are  nearly  equal 
throughout  the  Mach  number  regime.  These  trends 
are  true  for  the  other  flow  field  results. 


Finally,  Figure  IS  shows  volumetric 
efficiency,  V2/3/S  (where  V  is  the  total  volume  of 
the  waverider  and  S  is  its  planform  area)  versus 
Mach  number  for  the  optimum  waveriders.  This 
figure  is  shown  to  indicate  the  order  of  magnitude 
of  the  volumetric  efficiency  for  the  present  class 
of  waverider  vehicles.  In  general,  the  power-law 
body  derived  waveriders  appear  to  be  more  volume 
efficient  than  their  conically  derived  counter¬ 
parts,  but  the  trend  is  not  consistent. 

C,  Maximum  L/0  Versus  Minimum  Total  Drag 

This  section  compares  the  waveriders  that  are 
optimized  for  maximum  L/D  with  those  optimized  for 
minimum  total  drag.  Maximizing  for  L/0  can  be 
thought  of  as  searching  for  the  most  efficient 
hypersonic  cruising  configuration,  while  mini¬ 
mizing  Cg  corresponds  to  finding  the  best  acce¬ 
lerator  type  hypersonic  vehicle.  Figure  19 
compares  the  configurations  of  two  cone  flow 
derived  waveriders  at  Mach  6  for  these  two  types 
of  optimization.  The  same  comparison  is  made  for 
waveriders  designed  from  flow  over  a  3/4  power-law 
body,  at  a  Mach  number  of  14,  in  Figure  20.  It  is 
evident  tha*  optimizing  for  maximum  L/0  and  mini¬ 
mizing  for  total  drag  lead  to  quite  different 
waverider  configurations.  For  the  cone-flow  Mach 
6  case  shown  in  Figure  19,  the  value  of  L/0  is 
7.175  when  maximizing  L/0  and  7.187  when  mini¬ 
mizing  drag.  The  value  of  the  drag  coefficient, 
Cq,  for  the  drag  minimization  case  (Cg  =  0.004222) 
is  lower  than  for  the  case  of  maximizing  L/D 
(Cg  =  0.004694).  The  corresponding  numbers  for 
the  Mach  14,  power-law  body  case  are  L/0  equal  to 
6.484  and  Cg  of  0.003095  for  the  maximization  of 
L/0  and  L/0  equal  to  6.533  and  CO  of  0.003033  for 
the  total  drag  minimization.  These  are  the  usual 
trends  seen  in  the  data  for  all  of  the  cases  run. 

D.  Study  Concerning  the  Nature  of  the 

Boundary  Layer 

A  study  was  made  to  determine  the  effect  of 
using  an  all  laminar,  an  all  turbulent,  and  a 
transiting  boundary  layer  in  the  design  of  an 
optimized  waverider.  Figure  21  shows  represen¬ 
tative  results  for  optimization  using  flow  over  a 
1/2  power-law  body  at  a  Mach  number  of  6.  The 
figure  shows  perspective  views  of  the  resulting 
waveriders  when  the  different  types  of  boundary 
layers  are  assumed.  The  value  of  L/0  is  9.97  for 
the  all  lamina'-  case,  7.35  for  the  all  turbulent 
case,  and  7.46  for  the  transisting  case.  The  low 
value  of  L/0  for  the  transisting  flow  suggests 
that  much  of  the  waverider  surface  is  covered  with 
turbulent  flow.  The  data  for  this  case  confirms 
this  conclusion,  indicating  an  early  transition 
from  laminar  to  turbulent  flow. 

F.  Computing  Time 


IV.  CONCLUSIONS 

The  following  conclusions  are  made  from  the 
present  waverider  optimization  study: 

1.  The  present  results  are  the  first  examples  of 
viscous  optimizated  waveriders  designed  from 
axisymmetric  flow  fields,  such  as  the  hyper¬ 
sonic  flow  fields  over  minimum  drag  bodies. 

2.  The  waveriders  designed  from  flow  fields  over 
power- law  bodies  show  promise,  in  terms  of 
high  L/D  and  minimum  drag,  when  compared  with 
conical-flow  waveriders  or  other  hypersonic 
configurations. 

3.  Optimization  for  a  hypersonic  cruise  aircraft 
(maximizing  L/D)  results  in  a  distinctly  dif¬ 
ferent  configuration  than  optimization  for  a 
hypersonic  accelerator  aircraft  (minimizing 

Cg)* 

4.  The  transition  model  used  in  the  present  study 
results  in  turbulent  flow  over  much  of  the 
waverider  surface.  The  nature  of  the  boundary 
layer  has  a  large  effect  on  the  aerodynamic 
characteristics  of  the  vehicle.  Better  pre¬ 
diction  methods  for  transition  in  hypersonic 
flows  is  badly  needed. 

V.  ACKNOWLEDGEMENT 

This  work  has  been  supported  by  the  Minta 
Martin  Aeronautical  Fund  of  the  University  of 
Maryland.  Also  in  support  of  this  work,  the 
authors  have  enjoyed  a  close  working  relationship 
with  the  McDonnell  Douglas  Aircraft  Company  and 
the  Supersonic/Hypersonic  Branch  at  the  NASA 
Langley  Research  Center. 

REFERENCES 

1.  Anderson,  John  D.,  Jr.,  "A  Survey  of  Modern 
Research  in  Hypersonics,"  A1AA  Paper  84-1578, 
invited  paper  given  at  the  AIAA  17th  Fluid 
Dynamics,  Plasma  Dynamics,  and  Lasers 
Conference,  Snowmass,  Colorado,  June  25-27, 
1984. 

2.  Bowcutt,  Kevin  G. ,  Anderson,  John  D.,  Jr.,  and 
Capriotti,  D.P.,  "Viscous  Optimized  Hypersonic 
Waveriders,"  AIAA  Paper  87-0272,  presented  at 
the  AIAA  25th  Aerospace  Sciences  Meeting, 

Reno,  Nevada,  January  12-15,  1987. 

3.  Bowcutt,  K.G.,  Anderson,  J.D.,  Jr.,  and 
Capriotti,  D.P.,  "Numerical  Optimization  of 
Conical  Flow  Waveriders  Including  Detailed 
Viscous  Effects,"  in  Aerodynamics  of 
Hypersonic  Lifting  VEhicles,  AGARD-CPP-428, 
March  1987,  pp.  27/1-27/23. 


As  stated  earlier,  100  optimization  levels  or  4.  Wallace  Sawyer  and  Patrick  Johnston,  NASA 
iterations  are  used  in  the  optimization  of  a  Langley  Research  Center,  private  com- 

waverider  configuration.  This  requires  approxi-  munication,  1985. 

mately  5  to  15  minutes  of  CPU  time  on  a  UNIVAC 

1190  computer,  depending  on  the  type  of  flow  field  5.  Bowcutt,  Kevin  G. ,  Optimization  of  Hypersonic 
generating  body  used.  The  short  running  time  of  Waveriders  Derived  from  Cone  Flows--lncluding 

the  waverider  optimization  program  makes  it  Viscous  Effects.  Ph.D.  Dissertation, 

possible  to  conduct  parametric  studies  with  these  Department  of  Aerospace  Engineering, 

class  of  vehicles.  University  of  Maryland,  College  Park, 

Maryland,  1986. 


6.  Corda,  Stephen,  Ph.D.  Dissertation,  Department 
of  Aerospace  Engineering,  University  of 
Maryland,  College  Park,  Maryland,  to  be 
completed  December  1987. 

7.  Nonweiler,  T.R.F.,  “Aerodynamic  Problems  of 
Manned  Space  Vehicles,"  Journal  of  the  Royal 
Aeronautical  Society,  Vol .  63,  1959, 

pp.  521-528. 

8.  Seddon,  J.  and  Spence,  A.,  "The  Use  of  Known 
flow  Fields  as  an  Approach  to  the  Design  of 
High  Speed  Aircraft,"  in  Hypersonic  Boundary 
Layers  and  Flow  Fields,  AGARD  Cl>  No.  30,  May 
1968,  pp.  10/1-10/21. 

9.  Towend,  L.H.,  "Research  and  Design  for  Lifting 
Reentry,"  Progress  in  Aerospace  Sciences, 

Vol.  18,  19757 "pp."  i-SO. 

10.  Jones,  0.5. ,  Moore,  K.C.,  Pike,  J.,  and  Roe, 
P.L.,  "A  Method  for  Designing  Lifting 
Confi gurations  for  High  Supersonic  Speeds, 
Using  Axisymmetric  Flow  Fields," 
Inqenieur-Archiv,  Vol.  37,  pp.  56-72. 

11.  Rasmussen,  M.L.,  "Waverrider  Configurations 
Derived  from  Inclined  Circular  and  Elliptic 
Cones,"  J.  of  Spacecraft  and  Rockets,  Vol.  17, 
No.  6,  Nov. -Dec.  1980,  pp.  537-545. 

12.  Rasmussen,  M.L.,  and  Clement,  L.W. , 
"Cone-Derived  Waveriders  with  Longitudinal 
Curvature,"  AIAA  Paper  No.  84-2100,  1984. 

13.  Cole,  J.O. ,  and  Zien,  T.F. ,  “A  Class  of 
Three-Dimensional  Optimum  Hypersonic  Wings," 
AIAA  Journal,  Vol.  7,  No.  2,  Feb.  1969, 

pp.  264-271. 

14.  Kim,  B.S.,  Rasmussen,  M.L.,  and  Jischke,  M.C., 
"Otimization  of  Waverider  Configurations 
Generated  from  Axisymmetric  Conical  Flows," 
AIAA  Paper  32-1299,  1982. 

15.  Nelder,  J.A.,  Mead,  R. ,  "A  Simplex  Method  for 
Function  Minimzation",  Computer  Journal , 

Vol.  7,  Jan.  1965,  pp.  308-313. 

16.  Zucrow,  M.J.,  and  Hoffman,  J.O.,  Gas  Dynamics, 
John  Wiley  and  Sons,  New  York,  1972. 

17.  Anderson,  John  D.,  Jr.,  Modern  Compressible 
Flow  with  Historical  Perspective,  McGraw-Hill 
Book  co..  New  York,  1982. 

18.  Eckert,  E.R.G.,  "Engineering  Relations  for 
Heat  Transer  and  Friction  in  High-Velocity 
Laminar  and  Turbulent  Boundary-Layer  Flow  Over 
Surfaces  With  Constant  Pressure  and 
Temperature,"  Transactions  of  the  ASME, 

Vol.  78,  No.  6,  August  1956. 

19.  DiCristina,  V.,  “Three-Dimensional  Laminar 
Boundary  Transition  on  a  Sharp  8  Degree  Cone 
at  Mach  10",  AIAA  Journal ,  Vol.  8,  No.  5,  May 
1970,  p.  855. 

20.  Van  Wie,  O.M.,  An  Application  of  Computational 
Fluid  Dynamics  to  the  Design  of  Optimum  Ramjet 
and  Powered  Missile  Components,  Pn.O. 
Dissertation,  Department  of  Aerospace 


Engineering,  University  of  Maryland,  College 
Park,  Maryland,  1985. 

21.  Fetterman,  D.E.,  "Favorable  Interference 
Effects  on  Maximum  Lift-Drag  Ratios  on 
Half-Cone  Delta-Wing  Configurations  at  Mach 
6.86,"  NASA  TN  D-2942,  August  1965. 


TABLE  1 

Experimental  and  Numerical  Calculations 
of  tne  Hypersonic  Flow  over  a  Cone-Wing 


H,  «  6, 

,86,  ReL  ■  1. 

■  4J.106,  L  *  6, 

,404  in. 

a  *  4° 

9cone  • 

‘delta  *  «*' 

,  T,  -  59.2  X, 

T0  •  617 

K 

Case 

'.ill1*1 

CL 

CD 

I/O 

experiment 

? 

.012 

.0045 

2.67 

numerical 
(freest ream  upper 
surface) 

59.2 

617. 

.0126 

.0126 

.00329 

.00294 

3.83 

4.29 

numerical 
(1*  compression 
upper  surface) 

59.2 

617. 

.0126 

.0126 

.00198 

.00115 

1.05 

1.10 

Moch  Number 


FIG.  1:  Maximum  lift-to-drag  ratio  comparison 
for  various  hypersonic  configurations 
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FIG.  2:  Nonweiler  or  caret  wing  (from  Ref.  8) 
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FIG.  3:  Cone-flow  wing  (from  Ref.  8) 
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ABSTRACT 


A  family  of  optimized  hypersonic  waveriders 
is  generated  and  studied  wherein  detailed  viscous 
effects  are  included  within  the  optimization  pro¬ 
cess  itself.  This  is  in  contrast  to  previous 
optimized  waverider  work,  wherein  purely  inviscid 
flow  is  used  to  obtain  the  waverider  shapes.  For 
the  present  waveriders,  the  undersurface  is  a 
streamsurface  of  an  inviscid  conical  flowfield, 
the  upper  surface  is  a  streamsurface  of  the 
inviscid  flow  over  a  tapered  cylinder  (calculated 
by  the  axisymmetric  method  of  characteristics), 
and  the  viscous  effects  are  treated  by  integral 
solutions  of  the  boundary  layer  equations. 
Transition  from  laminar  to  turbulent  flow  is 
Included  within  the  viscous  calculations.  The 
optimization  is  carried  out  using  a  non-linear 
simplex  method.  The  resulting  family  of  viscous 
hypersonic  waveriders  yields  predicted  high 
values  of  lift/drag,  high  enough  to  break  the 
"L/0  barrier"  based  on  experience  with  other 
hypersonic  configurations.  Moreover,  the 
numerical  optimization  process  for  the  viscous 
waveriders  results  in  distinctly  different 
shapes  compared  to  previous  work  with  inviscid 
-designed  waveriders.  Also,  the  fine  details 
of  the  viscous  solution,  such  as  how  the 
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shear  stress  is  distributed  over  the  surface,  and 
the  location  of  transition,  are  crucial  to  the 
details  of  the  resulting  waverider  geometry. 

I.  INTRODUCTION 

Over  the  past  few  years.  Interest  in  all 
aspects  of  hypersonic  flight  has  grown  explosi¬ 
vely,  driven  by  new  vehicle  concepts  such  as  the 
National  Aerospace  Plane  (NASP),  aero-assisted 
orbital  transfer  vehicles  (AOTV's),  the  hyper¬ 
sonic  transport  (the  "Orient  Express"),  and 
hypersonic  missiles,  to  name  just  a  few.  An 
extended  discussion  of  these  concepts,  as  well  as 
a  survey  of  hypersonic  aerodynamic  research 
contrasting  the  "old"  with  the  "new"  hypersonics, 
is  given  in  Ref.l.  Hence  no  further  elaboration 
will  be  given  here. 

The  present  paper  deals  with  a  class  of 
advanced  hypersonic  lifting  configurations.  To 
help  understand  the  motivation  for  the  present 
work,  the  following  background  is  given.  For  a 
lifting  aerodynamic  body,  it  is  well-known  that 
high  maximum  lift-to-drag  ratios,  (L/D)max,  are 
very  difficult  to  obtain  at  hypersonic  speeds, 
due  to  the  presence  of  strong  shock  waves  (hence 
high  wave  drag)  and  massive  viscous  effects.  At 
supersonic  and  hypersonic  speeds,  the  most  effi¬ 
cient  lifting  surface  is  the  infinitely  thin  flat 
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plate;  the  Inviscld  hypersonic  aerodynamic  pro¬ 
perties  of  a  flat  plate  are  shown  as  the  solid 
curves  In  Fig.l,  based  on  the  Newtonian  limit  of 
free  stream  Mach  number  H»-*»  and  y  =  Cp/Cy-*,1. 
Note  that  L/D  theoretically  approaches  infinity 
as  the  angle-of-attack,  a,  approaches  zero.  In 
reality,  viscous  effects  will  cause  L/0  to  peak 
at  low  values  of  a,  and  to  go  to  zero  as  a->0. 

This  is  illustrated  by  the  dashed  line  in  Fig.l, 
which  shows  the  variation  of  L/0  modified  by  skin 
friction  as  predicted  by  a  reference  temperature 
method  described  in  Ref.  2.  Although  the  infinit- 
ly  thin  flat  plate  shown  in  Fig.l  is  the  most 
effective  lifting  surface  aerodynamically,  it  is 
the  least  effective  in  terms  of  volume  capacity. 

It  goes  without  saying  that  all  practical  flight 
vehicles  must  have  a  finite  volume  to  carry  fuel, 
payload,  etc.  Hence,  the  flat  plate  results,  al¬ 
though  instructive,  are  of  academic  interest 
only.  In  contrast.  Fig.  2  shows  values  of 
(L/0),fl|X  versus  the  volume  parameter  Vz/1 2 3/S  for 
several  generic  hypersonic  configurations,  obtain- 
from  Ref.  3.  Here,  V  is  the  body  volume  and  S  is 
the  planform  area.  Note  from  Fig. 2  that  typical 
hypersonic  values  of  ( L / D ) rna x  range  from  4  to  6 
for  such  lifting  bodies  at  the  conditions  shown. 
These  values  are  also  typical  of  the  hypersonic 
transport  configuration  studied  in  Ref.  4.  Clear¬ 
ly,  values  Of  (L/0)max  for  hypersonic  vehicles 
are  substantially  lower  than  those  for  conven- 
al  subsonic  and  low  supersonic  airplanes.  (For 
example,  the  maximum  L/0  values  for  the  World  War 
II  Boeing  B-29  and  the  contemporary  General  Dy¬ 
namics  F-lll  are  16.8  and  15.8  respectively,  as 
obtained  from  Ref.  5).  Indeed,  as  increases 
across  the  supersonic  and  hypersonic  regimes, 
there  is  a  general  empirical  correlation  for 
( L/D )max  ba$ed  on  actuaWl  ight  vehicle  experi¬ 
ence,  given  by  Kuchemann  : 


4  (MO) 

<L'D>max  *  -* - 

m 

This  variation  is  shown  as  the  solid  curve  in 
Fig. 3.  41so  shown  are  a  number  of  data  points 
for  various  previous  hypersonic  vehicle  configu¬ 
rations  at  various  Reynolds  numbers  (the  open  sym¬ 
bols),  as  well  as  new  results  from  the  present  in¬ 
vestigation  (the  solid  symbols).  Fig.  3  is  pivo¬ 
tal  to  the  present  paper,  and  will  be 
discussed  at  length  in  subsequent  sections.  How¬ 
ever,  at  this  stage  in  our  discussion.  Fig. 3  is 
used  to  illustrate  only  the  following  aspects: 

1.  The  solid  curve  represents  a  type  of  "L/D 
barrier"  for  conventional  hypersonic  vehi¬ 
cles,  which  is  difficult  to  break. 

2.  Data  for  conventional  hypersonic  vehicles, 
shown  as  the  open  circles, 

form  an  almost  random  "shotgun"  pattern 
which,  for  the  most  part,  falls  below  the 
solid  curve. 

(The  numbers  adjacent  to  these  open 
circles  pertain  to  specific  reference  num¬ 
bers  itemized  in  Ref.  7,  which  should  be 
consulted  for  details.) 

3.  The  solid  symbols  pertain  to  the  present 
study,  and  represent  a  new  class  of  hyper¬ 
sonic  configurations  which  break  the  "L/D 
barrier."  These  configurations  are  conical 
flow  waveriders  that  are  optimized  with  de¬ 
tailed  viscous  effects  included  directly  in 
the  optimization  process. 


To  help  understand  the  contribution  made  by 
the  present  work,  let  us  briefly  review  the  gener¬ 
al  concept  of  waveriders.  In  1959,  the  design  of 
three-dimensional  hypersonic  vehicles  which  support 
planar  attajh|d  shock  waves  was  Introduced  by 
Nonwei Her,  ’  who  hypothesized  that  streamsurfaces 
from  the  flow  behind  a  planar  oblique  shock  could 
be  used  as  supersonic  lifting  surfaces.  This  led 
to  a  class  of  vehicles  with  a  caret-shaped  trans¬ 
verse  cross-section  and  a  delta  planform--  the  so- 
called  caret  wing  as  shown  in  Fig. 4.  Here,  the 
body  surface  is  generated  by  stream  surfaces  behind 
a  planar  oblique  shock  wave.  The  shock  wave  is 
attached  to  the  sharp  leading  edges  at  the  design 
Mach  number,  and  hence  no  flow  spillage  takes 
place  around  the  leading  edge.  The  lift  is  high 
due  to  the  high  pressures  behind  a  two-dimensional 
planar  shock  wave,  exerted  on  the  lower  surface  of 
the  vehicle.  Because  the  body  appears  to  be  rid¬ 
ing  on  top  of  the  attached  shock  wave,  it  is 
called  a  "waverider".  The  aerodynamic  advantages 
of  such  waveriders  are  listed  in  Ref.  I,  and  are 
discussed  in  great  detail  in  Refs.  7  and  10.  In 
short,  without  repeating  the  details  here,  at  a 
given  lift  coefficient,  caret  waveriders  theore¬ 
tically  operate  at  higher  L/D  values  than  other 
hypersonic  conf igurations . 

Expanding  on  this  philosophy,  other  types  of 
flowfields  can  be  used  to  generate  waveriders.  For 
example,  any  streamsurface  from  the  supersonic 
flow  over  an  axisymmetric  body  can  be  used  to 
generate  a  waverider  with  an  attached  shock  wave 
along  its  complete  leading  edge.  Work  on  such 
waveriders  was  first  carried  out  in  Britain,  as 
nicely  summarized  in  Ref.  11,  where  the  flow  over 
a  right-circular  cone  at  zero  degrees  angle  of 
attack  is  used  to  generate  a  class  of  "conical 
flow"  waveriders.  Still  later,  waveriders  were 
generated  from  inclined  circular  and  elliptic 
cones,  and  axisymmetric  bodies  with  longitudinal 
curvature  by  Rasmussen  et  al .  ’  ,  using  hyper¬ 

sonic  small  disturbance  theory.  This  work  was 
further  embelished  by  the  search  for  optimized 
waverider  shapes.  for  example  Cole  and  Zien 
found  optimized  waveriders  derived  from  axisym¬ 
metric  bodies  with  longitudinal  curvature  by  using 
hypersonic  small  disturbance  theory  to  generate 
inviscid  flow  solutions,  and  then  utilizing  the 
calculus  of  variations  to  obtain  the  optimum 
waverider  shapes.  Later,  Kim  et  al.  used  the  same 
philosophy  to  derive  optimum  waveriderSj from 
flowfields  about  unyawed  circular  cones  ,  and 
yawed  circular  and  elliptic  cones  .  In  Refs.  14- 
16,  the  advantage  of  hypersonic  small  disturbance 
theory  is  that  analytic  expressions  are  obtained 
for  surface  pressure  distributions,  hence  lift  and 
wave  drag,  thus  enabling  the  application  of  the 
calculus  of  variations  for  optimization. 

Unfortunately,  to  date  the  potential  superiority 
of  waveriders  as  hypersonic  high  L/D  shapes  has 
not  been  fully  demonstrated,  either  in  the  wind 
tunnel  or  in  flight.  A  basic  problem  arises 
because  of  the  tendency  for  waveriders  to  have 
large  wetted  surface  areas,  which  leads  to  large 
friction  drag.  All  previous  waverider  optimiza¬ 
tion  work  (such  as  Refs.  14-16)  has  been  based  on 
the  assumption  of  inviscid  flow,  after  which  an 
estimate  of  skin  friction  for  the  resulting  con¬ 
figuration  is  sometimes  added.  As  a  result,  the 
real  aerodynamic  performance  of  the  resulting  opti¬ 
mum  configuration  usually  falls  short  of  its  expec¬ 
tations. 
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The  purpose  of  the  present  work  is  to  remove 
this  deficiency.  In  particular,  a  series  of  con¬ 
ical-flow  generateo  waveriders  are  optimized  for 
maximum  l/D  wherein  detailed  viscous  effects  (in¬ 
cluding  boundary  layer  transition)  are  included 
within  the  optimization  process  itself.  This 
leads  to  a  new  class  of  waveriders  where  the  opti¬ 
mization  process  is  trying  to  reduce  the  wetted 
surface  area,  hence  reducing  skin  friction  drag, 
while  at  the  same  time  maximizing  L/D.  Because 
detailed  viscous  effects  can  not  be  couched  in 
simple  analytical  forms,  the  formal  optimization 
methods  based  on  the  calculus  of  variations  can 
not  be  used.  Instead,  in  the  present  work  a 
numerical  optimization  technique  is  us^d,  based  on 
the  simplex  method  by  Nelder  and  Mead1  .  By  using 
a  numerical  optimization  technique,  other  real 
configuration  aspects  can  be  included  in  the  ana¬ 
lysis  in  addition  to  viscous  effects,  such  as 
blunted  leading  edges,  and  an  expansion  upper  sur¬ 
face  (in  contrast  to  the  standard  assumption  of  a 
free  stream  upper  surface,  i.e.,  an  upper  surface 
with  all  generators  parallel  to  the  freestream 
direction).  The  results  of  the  present  study  lead 
to  a  new  class  of  waveriders,  namely  "viscous 
optimized"  waveriders.  Moreover,  these  waveriders 
appear  to  produce  relatively  high  values  of  (L/D), 
as  will  be  discussed  in  subsequent  sections. 

II.  ANALYSIS 

For  the  present  waverider  configurations,  the 
following  philosophy  is  followed: 

1.  The  lower  (compression)  surface  is  generated 
by  a  streamsurface  behind  a  conical  shock 
wave.  The  inviscid  conical  flowfield  is 
obtained  from  the  numerical  solution  of  the 
Taylor-Maccol 1  equation,  derived  for  example 
in  Ref.  18. 

2.  The  upper  surface  is  treated  as  an  expansion 
surface,  generated  in  a  similar  manner  from 
the  inviscid  flow  about  a  tapered,  axisym- 
metric  cylinder  at  zero  angle  of  attack,  and 
calculated  by  means  of  the  axisymnetric 
method  of  characteristics. 

3.  The  viscous  effects  are  calculated  by  means 
of  an  integral  .boundary  layer  analysis 
following  surface  streamlines,  including 
transition  from  laminar  to  turbulent  flow. 

4.  Blunt  leading  edges  are  included  to  the 
extent  of  determining  the  maximum  leading 
edge  radius  required  to  yield  acceptable 
leading  edge  sgfface  temperatures,  and  then 
estimating  the  leading  edge  drag  by  modified 
Newtonian  theory. 

5.  The  final  waverider  configuration,  optimized 
for  maximum  L/D  at  a  given  Mach  number  and 
Reynolds  number  with  body  fineness  ratio  as 

a  constraint,  is  obtained  from  the  numerical 
simplex  method  taking  into  account  all  the 
effects  itemized  in  1-4  above  within  the 
optimization  process  Itself. 

The  following  paragraphs  describe  each  of  the 
above  Items  In  more  depth;  for  a  highly  detailed 
discussion,  see  Ref.  7. 


A.  Inviscid  Flow  --  Lower  Surface 

The  waverider's  lower  surface  is  generated 
from  a  streamsurface  behind  a  conical  shock  wave 
supported  by  a  hypothetical  right  circular  cone 
at  zero  angle  of  attack.  The  hypothetical  cone 
and  its  flowfield  is  shown  in  Fig.  5,  where  0C  is 
the  cone  semi-angle  and  es  is  the  wave  angle.  The 
inviscid  conical^low  is  obtained  from  the  Taylor- 
Maccol  1  equation19 
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solved  by  a  standard  Runge-Kutta,  forth-order 
accurate  numerical  technique,  namely  the  RKF45 
algorithm  obtained  from  Ref.  19.  In  Eq.(l),  Vr  is 
the  component  of  nondimensional  flow  velocity 
along  a  conical  ray,  e  is  the  angle  of  the  ray 
referred  to  the  cone  axis,  and  y  is  the  ratio  of 
specific  heats. 

Any  streamsurface  from  this  flowfield  can 
represent  the  wing  undersurface  of  a  waverider, 
as  shown  in  Fig.  6.  (For  purposes  of  illustra¬ 
tion,  Fig.  6  also  shows  the  waverider  upper  sur¬ 
face  as  a  freestream  surface,  but  this  is  only 
one  of  many  possible  choices  for  the  upper 
surface.)  Any  particular  undersurface  is  uni¬ 
quely  defined  by  the  intersection  of  the  conical 
flow  streamsurface  with  the  conical  shock  wave, 
as  shown  by  the  curve  labeled  "leading  edge"  in 
Fig.  7.  Let  us  examine  Fig. 7  more  closely.  It 
is  a  front-view  of  the  hypothetical  conical 
flowfield,  illustrating  the  cone  apex  at  the 
center,  and  both  the  cone  base  and  shock  base  at 
some  arbitrary  distance  downstream  of  the  apex. 
Consider  a  curve  in  this  front-view,  lying  below 
the  apex  (or  even  including  the  apex),  as  shown 
by  the  curve  labeled  "leading  edge".  Now 
construe  this  curve  as  a  trace  on  the  conical 
shock  wave  itself,  and  visualize  streamlines 
trailing  downstream  from  this  trace;  the 
resulting  streamsurface  is  the  waverider  under¬ 
surface  sketched  in  Fig.  6.  Indeed,  the  curve 
labeled  "leading  edge”  in  Fig.  7  is  simply  the 
forward  projection  of  the  waverider  leading  edge 
on  the  cross-flow  (x-y)  plane.  This  curve  Is 
treated  as  completely  general,  except  for  the 
constraints  that  it  be  symmetric  about  the  y-z 
plane,  and  that  it  lie  entirely  below  the  x-z 
plane  to  ensure  that  the  waverider  undersurface 
is  a  compression  surface.  Also  in  Fig.  7,  note 
the  curve  labeled  "trailing  edge".  This  is  the 
intersection  of  the  particular  conical  flow 
streamsurface  with  the  plane  of  the  shock  base, 
and  it  represents  the  bottom  surface  trailing 
edge  of  the  waverider.  This  Is  the  shape  of  the 
bottom  of  the  waverider  base,  as  sketched  in  Fig. 
T Returning  to  Fig. 7,  the  area  between  the 
"leading  edge"  and  “trailing  edge"  curves  Is  the 
forward  projection  of  the  entire  waverider 
compression  surface.  Moreover,  the  dashed  line 
emanating  from  the  cone  apex  in  Fig.  7  is  the 
forward  projection  of  a  conical  flow  streamline; 
hence,  that  portion  of  the  dashed  line  contained 
between  the  ^leading  edge"  and  "trailing  edge" 
curves  Is  the  projection  of  a  particular 
streamline  along  the  waverider  undersurface,  from 
the  leading  edge  to  the  trailing  edge. 


,  B .  Inviscid  Flow  --  Upper  Surface 

« 

<  In  most  previous  waverider  work,  the  upper 

surface  is  treated  as  a  freestreara  surface,  as 
“  illustrated  in  Fig.  6.  Here,  the  upper  surface 

pressure  is  freestream  pressure,  p..  However,  if 
i  the  upper  surface  is  made  an  expansion  surface, 

i  where  p<p„,  then  a  small  but  meaningful  contribu- 

|  tion  to  I/O  can  be  obtained.  This  approach  is 

.  taken  here.  Similar  to  the  philosophy  for  the 

.  lower  surface,  the  upper  surface  is  a  streamsur- 

face  “carved"  from  a  known  expansion  flow.  The 
hypothetical  expansion  body  chosen  here  is  a  cir¬ 
cular  cylinder  of  given  radius;  the  cyclinder  is 
|  aligned  parallel  to  the  flow  and,  at  some  point, 

!  is  tapered  parabol ical ly  to  a  smaller  radius.  The 

'  result  is  an  axisymmetric  expansion  flow,  where 

the  domain  of  expansion  is  bounded  by  a  freestream 
’  Mach  cone  centered  on  the  cy  linder  axis,  as  shown 

;  in  Fig.  B.  Parabolic  taper  was  chosen  because  it 

;  is  relatively  simple,  and  the  resulting  expansion 

body  slope  is  everywhere  continuous.  Once  the 
expansion  body  is  chosen,  it  remains  only  to 
geometrical ly  position  the  expansion  region  rela¬ 
tive  to  the  lower  surface,  choose  the  initial  and 
final  cylinder  radii,  solve  the  inviscid  expansion 
flow,  then  cut  a  st reamsurf ace  from  that  flow  to 
serve  as  the  waverider  upper  surface.  This  basic 
idea  was  first  -levelled  for  two-dimensional 
expansions  by  Flowej^  ,  and  later  for  axisymmetric 
expansions  by  Moore 

The  axisymmetric  flow  is  calculated  from  the 
axisynrietric  method  of  characteristics,  using  the 
two-st^p  predictor-corrector  iteration  scheme  of 
Ferri  .  The  details  involving  the  matching  of 
the  resulting  expansion  surface  with  the  conical 
flow  compression  surface  are  straight  forward,  but 
lengthy.  Considering  that  the  expansion  surface 
contributes  only  about  10%  to  the  value  of  (I/O), 
no  further  space  for  its  discussion  is  justified 
in  the  present  paper;  for  the  complete 
discussion,  see  Ref.  7. 

C.  Leading-Edge  Bluntness 

Waveriders,  by  design,  have  sharp  leading 
edges  that  support  attached  shock  waves. 

However,  for  flight  Mach  numbers  above  five,  the 
temperatures  for  sharp  leading  edges  will  exceed 
the  practical  limits  of  most  structural 
materials.  This  leads  to  the  need  for  blunt 
leading  edges  with  sufficiently  large  radii  such 
that  the  aerodynamic  heat  flux  is  reduced  to 
reasonable  levels.  However,  at  the  same  time  the 
leading  edge  radius  should  be  as  small  as 
possible  to  reduce  the  nose  drag. 

To  re^yce  the  required  leading  edge  radius, 
Nonwefler23  has  proposed  adding  conducting 
material  aft  of  the  leading  edge  to  transport 
thermal  energy  away  from  the  region  of  high  con¬ 
vective  heating  near  the  stagnation  or  attachment 
line,  and  conduct  it  downstream  to  areas  where 
convective  heating  is  lower,  and  excess  energy  can 
be  radiated  away  from  the  body.  Nonweiler  labeled 
this  theoretical  concept  as  a  "conducting  plate", 
which  is  somewhat  analogous  to  other  passive 
cooling  techniques,  such  as  heat  pipes.  Using 
Nonweiler's  basic  technique,  minimum  leading  edge 
radii  can  be  ascertained,  once  flight  Mach  number, 
freestream  conditions,  leading  edge  sweep, 
material  properties,  and  maximum  allowable  tem¬ 
perature  are  known. 


In  the  present  work,  Monweiler's  technique 
was  used  to  determine  the  leading  edge  radii  for 
waveriders  designed  for  Mach  nunbers  between  6 
and  25.  The  leading  edge  material  used  for  the 
calculations  was  ATJ  graphite,  chosen  because  it 
is  representative  of  materials  with  high  conduc 
tivity  and  high  melting  point  temperature.  Details 
of  this  technique  as  applied  to  the  present  work 
are  given  in  Ref.  7.  It  is  interesting  to  note 
that,  for  conditions  associated  with  the  typical 
flight  path  of  a  lifting  hypersonic  vehicle  enter¬ 
ing  the  earth's  atmosphere  at  Mach  25  and  dece¬ 
lerating  to  Mach  6  at  lower  altitude,  the  minimum 
leading  edge  diameters  ranged  from  6  to  2B.  mi  -- 
quite  small  in  comparison  to  a  typical  overall 
length  of,  say  50m.  Therefore,  the  present  waveri¬ 
ders  are  essentially  “aerodynamically  sharp”  from 
that  point  of  view.  Regardless  of  the  apparently 
small  amount  of  required  leading  edge  bluntness 
(from  the  aerodynamic  heating  point  of  view),  the 
present  waveriders  were  geometrically  altered  to 
accomodate  the  blunt  edge,  and  the  contribution  to 
aerodynamic  forces  on  the  waverider  were  estimated 
assuming  a  modified  Newtonian  pressure  distribution 
on  the  leading  edge. 

0 .  Viscous  Flow  Analysis 

A  major  aspect  of  the  present  investigation 
is  that  optimum  waverider  shapes  are  obtained 
wherein  detailed  viscous  effects  are  included 
within  the  optimization  process  itself.  These 
viscous  effects  are  calculated  by  means  of  two 
integral  boundary  layer  techniques,  described 
below.  In  all  cases,  the  boundary  layer  flow  is 
assumed  to  be  locally  two-dimensional,  following 
the  inviscid  upper  surface  and  lower  surface 
streamlines.  Both  laminar  and  turbulent  flow  are 
considered,  along  with  a  transition  region  based 
on  empirical  correlations. 

D . 1  Laminar  Analysis 

The  laminar  boundary  layer  calculations 
were  performed  using  Walz  1  Integral  method,  as 
described  in  Ref.  24.  The  method  requires  the 
solution  of  a  set  of  coupled  first-order  ordinary 
diffe  ential  equations  along  the  boundary  layer 
edge  streamlines.  These  equations  are  the  boun- 
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equations. 
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momentum  and  mechanical  energy 
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Note  that  in  the  above  equations,  primes  denote 
differentiation  with  respect  to  x,  here  repre¬ 
senting  the  boundary  layer  coordinate  in  the 
streamline  direction.  The  variables  in  Eqs. 

( 9) - ( 12)  are  defined  by  Walz  to  be 


H  *  «j/S2  =  bH12+r  Y  M2(W-6) 
a  =  1.7261  (W*-1.515)°-71:>8 

(4^U  Y-l  2  ~ 

b  =  =  1+r  Y  Me(W-8)(2-W) 


r  =  /Pr 


Taw(x)'Vx) 

Taw^e^ 


and  8  =  8yx 
where 

H12  *  4.0306  -  4.2845  (W*-l .515) 


ru 


0.3886 


Taw  =  Te  +  IT 


w*  =  ^sV^u  =  W/* 


*  =  1  + 


[hr™* 
iy  1,1 12~ 1 
~*r~ 


2-(«  )  /6  _  !-(«,)/«  . 

■  — B*2-  •  *  -W-  "••I 


0.8 


*  0.0144(2-W*)(2-8) 


(«1)u/«  *  0.420  -  ( W* - 1 .515) 


0.424W* 


0.555 


g  -  0. 324+0. 336(W*-1. 515) 


Bu  *  0.1564  ♦  2.1921  (W*-1.515) 


1.70 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


(25) 

(26) 

(27) 


1  l+r(Y)^I  1.16W*-1.072-e  (2H*-2.581)]j 


0.7 


l+r(2y^-)K^(  1-8  )] 


ZU-0.7 


(28) 


Careful  inspection  of  Eq.  (21)  reveals  that  it  is 
an  implicit  equation  for  W* ,  since  4>  is  a  complex 
function  of  W*.  Therefore,  in  practice,  a  numeri¬ 
cal  zero-finding  routine  is  used  on  Eq.  (21)  to 
find  the  value  of  W*  that  yields  the  known  value  of 
W.  Walz  suggested  an  approximation  that  would 
allow  closed  form  solution  of  Eq.  (21),  however, 
the  present  authors  have  found  that  it  made  boun¬ 
dary  layer  calculations  blow  up  when  used  for  a 
Mac h  six  flat  plate  test  case.  Hence,  the 
suggested  approximation  was  discarded  in  favor  of 
the  zero-finding  approach.  For  more  details  on  the 
numerical  solution  of  these  integral  boundary  layer 
equations,  see  Ref.  7. 

D.2  Turbulent  Analysis 

If  and  when  boundary  layer  transition  is 
predicted,  turbulent  boundary  layer  calculations 
are  performed  using  the  inner  variable  integral 
method  of  White  and  Christoph  .  in  practice,  the 
method  requires  solution  of  one  of  two  first-order 
ordinary  differential  equations  along  the  boundary 
layer  edge  streamlines,  depending  upon  the  value  of 
the  parameter  X/X__„,  where 
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In  Eq.  ( 
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and  Taw  is  defined  by  Eq.  (20),  except  that  now 
the  recovery  factor  is  the  turbulent  value. 
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According  to  Ref.  24,  if  A/k 
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appl ies ,  where 


where  x^e  and  xtj  are  the  distances  along  a 
streamline  from'the  leading  edge  to  the  beginning 
and  end  of  transition,  respectively,  and  (Rex)^ 
the  local  Reynolds  number  at  the  beginning  of  tran¬ 
sition  obtained  in  the  present  analysis  from  Eq44  5j, 


The  variation  of  local  skin  friction  coef¬ 
ficient  within  the  transitional  region  (between  xte 
and  xtj)  is  assumed  to  be  a  linear  combination  of 
the  laminar  (cp^)  and  turbulent  (cfj)  values  that 

would  have  existed  if  the  boundary  layer  were 
completely  laminar  or  turbulent,  respectively.  The 
transitional  friction  coefficient,  CfTr<,  is  thus 
related  to  Cf^  and  Cf-j.  oy: 


CfTR  =  (1'°  +  5  C'T  (4?) 

where  $  is  a  weighting  factor  (a  function  of  x) 
inspired  by  Emmons  (as  discussed  in  Ref. 29).  For 
the  present  investigation,  the  following  expres¬ 
sion  for  5  is,  as  derived  in  detail  in  Ref.  7: 


=  ]_p-3texp[^-;-(Rex)?;?(x-xtj)]-ll2  ^ 


:(*)  1  l-e 


f*  =  (2.434a  r  1.4432'  )exp(  -44  z*1 )  (41) 

g*  *  1  -  2.32  +  1.76z3  (42) 

z=l-  X/X  (42) 

and  the  primes  ienote,  as  in  the  laminar  case, 
di f ferent i at i on  with  respect  to  the  streamline 
coorjiojte,  x.  For  more  details  concerning  the 
numerical  solution  of  these  equations,  see  again 
Ref.  7. 


It  is  not  possible  within  the  current  state-of- 
the-art  to  evaluate  the  accuracy  of  these  tran¬ 
sition  correlations.  After  a  study  of  the  existing 
literature,  the  present  authors  feel  that  the  above 
relations  form  a  practical  method  for  simulating 
transition  within  the  goals  of  the  present  study. 
They  provide  a  mechanism  for  assessing  the  effect 
of  transition  on  optimum  waverider  shapes;  indeed, 
as  discussed  in  the  results,  one  series  of  numeri¬ 
cal  experiments  is  conducted  wherein  the  transition 
location  is  varied  as  a  parameter. 


0 .  3  Transition  Analysis 


E .  Aerodynamic  Forces 


The  prediction  of  transition  from  lami¬ 
nar  to  turbulent  flow  at  hypersonic  speeds  is  a 
state-of-the-art  research  topic.  In  the  present 
analysis,  the  correlation  used  for  predicting  the 
onset  of  transition  is  based  on  two  sets  of  data: 
(1)  data  for  sharp  cones  at  zero  angle-of- 
attack  ;  and  (2)  data  for  winc^  with  blunt, 
swept  supersonic  leading  edges7  .  The  correla¬ 
tion  gives  local  transition  Reynolds  number  ReXt, 
as  a  function  of  local  edge  Mach  number,  Me,  asc 
fol lows7b : 

log.JRe.  )  =  6.421  exp  ( 1  .ROgxlO'V  •6'* 1 )  (44) 

* '  xt  e 


The  lift,  drag,  and  hence  L/D  is  calculated 
from  a  detailed  integration  of  the  local  surface 
pressures  and  shear  stress  over  the  waverider 
surface.  Consistent  with  wind  tunnel  practice 
as  well  as  other  literature,  base  drag  is  not 
included  in  the  present  results.  (For  example, 
all  the  data  shown  id  Fig.  3  do  not  include 
base  drag.)  This  is  done  to  enable  a  rational 
comparison  with  other  data.  Moreover,  at  very 
high  Mach  number,  the  base  drag  becomes  a  small 
quantity  in  comparison  to  forebody  drag.  Details 
on  the  pressure  and  shear  stress  numerical 
integration  can  be  found  in  Ref.  7. 


In  turn,  this  value  of  transition  Reynolds  number 
is  modijfjed  for  wing  leading-edge  sweep,  as 
follows 

""V, 

(«•».) 


t  A=0 


0.787cos4,346A-0.7221e‘°*099U+0.9464 

(45) 


where  A  is  the  sweep  angle,  and  (Re„t)A=0  is 
obtained  from  Eq. (44). 

Once  the  onset  of  transition  has  been  pre¬ 
dicted,  the  extent  of  the  transition  region, 
hence  the  end  of  transition,  is  predicted  using  a 
relationship  developed  by  Harris  and  Blanchard7*, 
as  follows:  „  , 

xte  *  *tlll+5<Re*>;?’  i  <46> 


F.  Waverider  Optimization 

Once  a  specific  shape  for  the  forward  leading 
edge  projection  of  a  waverider  is  chosen,  (such  as 
shown  in  Fig.  7),  the  techniques  outlined  in  the 
previous  sections  can  be  used  to  generate  the  cor¬ 
responding  waverider  and  evaluate  its  lift-to-drag 
ratio  (I/O).  Finding  the  leading  edge  shape  that 
maximizes  L/D,  with  all  other  parameters  held 
fixed,  then  requires  an  optimization  scheme  that 
can  systematical ly  change  the  projected  leading 
edge  shape  in  search  of  the  one  that  yields  maximum 
l/D.  Unfortunately,  most  existing  optimization 
schemes  require  that  the  function  of  interest  have 
an  analytical  description  —  a  requirement  not  pos¬ 
sible  in  the  present  work.  There  is  one  scheme, 
however,  a  non-linear  simplex  mp^hod  for  function 
minimization  by  Nelder  and  Mead  ,  that  requires 
nothing  more  than  the  ability  to  numerically 
evaluate  the 


function.  This  scheme  has  been  used  in  the  pre¬ 
sent  work  to  find  optimum  waveriders. 

In  general,  the  scheme  of  Ref. 17  minimizes  a 
function  of  n  variables  by  comparing  values  of 
the  function  at  (n+1)  vertices  of  a  "simplex", 
then  replacing  the  vertex  with  the  highest  func¬ 
tion  value  by  another  point  determined  via  the 
logic  of  the  scheme.  As  a  result  of  the 
algorithm  logic,  "tne  simplex  adapts  itself  to 
the  local  landscape  [of  the  function  surface], 
elongating  down  long  inclined  planes,  changing 
direction  on  encountering  a  valley  at  an  angle, 
and  contracting  in  the  neighborhood  of  a 
minimum",  according  to  Ref. 17.  In  this  scheme, 
three  operations  --  reflection,  contraction  and 
expansion  —  are  used  to  modify  the  current 
simplex  in  an  attempt  to  replace  the  vertex 
having  the  highest  function  value  with  one  having 
a  lower  value.  Each  of  three  operations  replace 
one  or  more  o*  the  (n+1)  points  (Po.Rj ... . ,Pn) 
that  define  the  current  simplex  in  n-dimensional 
space  vi*  '  points  that  yield  progressively 
sm^' .er  fane. i  on  values  (fp.f^ ,. . . ,fn)  At  the  new 
vertex  points.  A  graphic  illustration  of  how  the 
method  works  is  shown  in  Fig.  9  for  a  hypotheti¬ 
cal  function,  f,  of  two  variables,  Cj  and  Cj.  In 
the  figure,  a  triangle  with  vertices  on  the  func¬ 
tion  sjrface  represents  a  possible  simplex.  In 
the  optimization  process,  the  triangle  (simplex) 
flip-flops  down  the  function  valley,  expanding  if 
possible  to  speed  up  the  process,  then  contract¬ 
ing  when  it  straddles  the  minimum. 

To  use  the  simplex  method  for  optimizing 
waverider  l/D,  the  shape  of  the  forward  projec¬ 
tion  of  the  leading  edge  must  be  parameteri zed  in 
some  general  way.  In  the  present  work,  five 
points  in  the  x-y  plane,  lying  inside  of  the 
shock  domain,  were  chosen  to  represent  the  for¬ 
ward  leading  edge  projection.  A  cubic  spline-fit 
through  the  five  points  is  then  used  to  generate 
a  continuous  leading  edge.  One  of  the  five 
points,  the  symmetry  plane  point,  is  constrained 
to  lie  on  the  y-axis,  hence  its  x-  value  is 
always  zero.  This  leaves  nine  variables,  the 
remaining  x  and  y  values  of  the  leading  edge  pro¬ 
jection  points,  for  the  optimization  routine  to 
manipulate  in  search  of  an  optimum  waverider.  A 
set  of  leading  edge  coordinates  thus  represent  a 
single  vertex  point, 

Pi  =  <x2-*3’*4*W-v2-y3**4‘y5)i  (49! 

of  the  required  simplex,  where  xj  *  0  as 
explained,  and  the  function  to  be  minimized  is 
the  negative  of  the  lift-to-drag  ratio 

VV  =  <-L/D’i  (50) 

Note  that  the  five  leading  edge  points  are  used 
to  define  only  half  of  the  projected  leading  edge 
shape,  since  the  other  half  is  constrained  by 
vehicle  symmetry  to  be  the  mirror  image  of  the 
first  half. 

With  nine  variables  (n=9),  ten  points  (hence 
ten  leading  edge  shapes)  must  be  chosen  to  create 
the  initial  simplex.  In  the  present  work,  five 
polynomials  of  the  form 

yte  *  C1  +  C2xte  +  C3X 


and  five  cosine  curves  of  the  form 


Jie 


(52) 


*  te  \ 
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were  used  to  describe  the  initial  leading  edge 
shapes;  the  constants  Ci,C2»...C6  being  varied  to 
generate  a  set  of  distinct  shapes.  An  example  of  a 
set  of  initial  leading  edge  shapes  is  shown  in  Fig. 
10  --  the  bold  line  representing  the  final  shape 
associated  with  the  optimum  waverider  for  this 
case.  Also  note  that  in  the  present  work,  100 
steps  of  the  optimization  routine  were  executed  for 
all  cases  run,  though  a  convergence  criterion  could 
have  been  implemented  as  described  in  Ref.  17.  It 
was  found  that  one-hundred  steps  provided  adequate 
convergence  for  engineering  accuracy  (*  10~3 -10"1* ) 
without  using  excessive  computer  resources  to 
generate  an  optimized  waverider. 


For  more  details  on  the  optimization  scheme, 
see  Ref.  7. 

III.  RESULTS  AND  DISCUSSION 

The  present  results  are  divided  into  four  sec¬ 
tions,  as  follows:  (1)  a  presentation  of  optimum 
waverider  shapes  and  aerodynamic  characteristics  at 
Mach  6  and  25,  representing  two  extremes  of  the 
hypersonic  flight  spectrum;  (2)  a  numerical  experi¬ 
ment  to  assess  the  impact  of  boundary  layer  tran¬ 
sition  on  the  optimized  waverider  shapes;  (3)  an 
assessment  of  the  need  to  account  for  detailed  sur¬ 
face  variations  of  shear  stress  in  contrast  to  the 
use  of  an  average  skin  friction  coefficient  during 
the  optimization  process;  (4)  an  examination  of  the 
question:  if  the  skin  friction  is  deleted  from  the 
present  analysis,  what  type  of  optimized  inviscid 
waverider  configuration  is  produced? 


Due  to  the  specialized  nature  of  any  waverider 
generation  analysis.  Including  the  present  one,  it 
is  difficult  to  obtain  a  direct  benchmark  com¬ 
parison  with  existing  data  in  order  to  verify  the 
integrity  of  the  current  results.  However,  with 
the  present  analysis,  it  is  possible  to  calculate 
the  aerodynamic  properties  of  a  half-cone  with  a 
flat  delta  wing  mounted  on  top;  in  this  case  the 
wing  will  have  a  sweep  angle  corresponding  to  the 
shock  angle  of  the  cone,  and  the  body  will  be  at 
zero  degrees  angle  of  attack.  This  specialized 
case  was  calculated  at  Mach  6.8  for  a  half-cone  of 
ec  =  3.67°,  and  the  corresponding  wing  sweep  angle 
of  81°.  The  result  is  given  as  the  flagged  solid 
square  in  Fig. 3.  This  is  to  be  partly  compared 
with  the  point  labeled  P2a,  which  was  obtained  from 
Ref.  3,  and  which  corresponds  to  a  similar  flat- 
top  half-cone,  delta  wing  model,  but  at  conditions 
of  maximum  l/D,  hence  at  some  positive  angle  of 
attack .  About  the  only  point  to  be  made  here  is 
that  the  calculated  L/D  at  zero  angle  of  attack  is 
lower  than  the  measured  (L/0)max  at  some  angle-of- 
attack  —  a  proper  qualitativeresult.  The 
measured  L/D  at  zero  angle-of-attack  is  not  pre¬ 
sented  in  Ref.  3;  however,  through  a  personal 
inquiry  to  Patrick  Johnston  at  NASA  Langley,  the 
present  authors  have  been  told  that  the  measured 
L/D  at  zero-angle-of-attack  was  2.7  —  about  eight 
percent  higher  than  the  value  of  2.5  calculated 
with  the  present  analysis.  This  is  a  reasonable 
comparison,  and  if  anything,  seems  to  Indicate  that 
the  present  aerodynamic  analysis  is  conservative. 
(Please  note  that  the  comparisons  discussed  above 
are  for  a  given  configuration,  not  an  optimized 
waverider;  hence  any  degree  of 
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of  validation  here  pertains  to  the  aerodynamic 
portion  of  the  analysis  and  not  to  the  present 
optimization  process  itself.) 

A .  Representative  Waveriders 

In  Ref.  7,  a  series  of  optimized  waveriders  is 
generated,  including  cases  at  M„  =  4,6,10,15,20 
and  25.  The  conditions  correspond  to  altitude- 
velocity  points  along  a  typical  entry  flight  tra¬ 
jectory  of  a  lifting  hypersonic  vehicle,  such  as 
an  aerospace  plane.  In  the  present  section,  only 
the  results  at  =  6  and  25  are  presented  as 
represeotat i ve  of  the  two  extremes  of  the  flight 
spectrum.  Ref.  7  shoul J  be  consulted  for  addi¬ 
tional  results. 

Fig.  11  gives  values  of  (L/0),  Cj_,  and  volu¬ 
metric  efficiency,  n  *  l/2/3/Sp,  for  waveriders 
optimized  at  different  assumed  wave  angles  for  the 
conical  shock.  To  unierstand  this  more  fully, 
consider  tne  conical  flow  field  associated  with  a 
given  conical  shock  wave,  say  9S  =  ll3.  For  this 
value  of  2$  an  optimum  waverider  shape  is  obtained 
[refer  agam  to  the  bold  curve  in  Fig.  10).  The 
resulting  character st ics  of  this  optimized  wave- 
rider  are  then  plotted  on  Fig.  11  for  *  11°. 
Tnis  process  is  repeated  for  other  values  of  95, 
say  12°,  13°,  and  14°.  For  each  value  of  9S,  an 
optimized  waverider  is  obtained,  and  its  charac¬ 
teristics  plotted  in  Fig.  11  as  the  open  symbols. 
(The  solid  symbols  will  be  discussed  later.) 

Hence,  Fig.  11  pertains  to  an  entire  series  of 
optimized  waveriders.  However,  note  that  the  (L/0) 
curve  itself  nas  a  maximum  (in  this  case  for  9^  = 
12°).  This  yields  an  "optimum  of  the  optimums", 
and  defines  tne  final  viscous  optimized  waverider 
at  M»  *  6  for  the  flight  conditions  shown  in  Fig. 
11.  The  front  views  of  the  optimum  shapes  at  each 
value  of  9S  are  shown  in  Fig.  12,  and  the  corres¬ 
ponding  perspective  views  are  shown  in  Fig.  13. 
Finally,  a  summary  three-view  of  the  best  optimum 
(the  “optimum  of  the  optimum")  waverider,  which 
here  corresponds  to  9S  *  12°,  is  given  in  Fig,  14. 
Also  in  Figs.  12-14,  the  lines  on  the  upper  and 
lower  surfaces  of  the  waveriders  are  inviscid 
streamlines.  Note  in  these  figures  that  the  shape 
of  the  optimum  waverider  changes  considerably  with 
es.  Moreover,  examining  (for  example)  Figure  14, 
note  the  rather  complex  curvature  of  the  leading 
edge  in  both  the  planform  and  front  views;  the 
optimization  program  is  shaping  the  waverider  to 
adjust  both  wave  drag  and  skin  friction  drag  so 
that  the  overall  L/D  is  a  maximum.  Indeed,  it  was 
observed  in  all  of  the  present  results  that  the 
best  optimum  shape  at  any  given  M  results  in  the 
magnitudes  of  wave  drag  and  skin  Triction  drag 
being  approximately  the  same,  never  differing  by 
more  than  a  factor  of  two.  For  conical  shock 
angles  below  the  best  optimum  {for  example  9c  = 

11°  in  Figs.  12  and  13),  skin  friction  drag  is 
greater  than  wave  drag;  in  contrast,  for  conical 
shock  angles  above  the  best  optimum  (for  example 
0S  *  13°  and  14°  in  Figures  12  and  13),  skin  fric¬ 
tion  drag  is  less  than  wave  drag.  (Note:  For  a 
hypersonic  flat  plate,  using  Newtonian  theory  and 
an  average  skin  friction  coefficient  ,  it  can 
readily  be  shown  that  at  maximumum  I/O,  the  wave 
drag  is  twice  the  friction  drag.) 

The  results  in  Figs.  11-14  pertain  to  M^  *  6. 
An  analogous  set  of  results  for  the  other  extreme 
of  the  lifting  hypersonic  flight  spectrum  at  M»  = 


25  is  given  in  Figs.  15-13.  The  aerodynamic  char¬ 
acteristics  of  optimum  waveriders  for  »  7°, 8°, 3° 
and  10°  are  given  as  the  open  symbols  Fig.  15  (the 
solid  symbols  will  be  discussed  later.)  The 
respective  front  views  are  shown  in  Fig.  16,  and 
perspective  views  in  Fig.  17.  Finally,  the  best 
optimum  Mach  25  waverider  (which  occurs  at  0S  = 

93)  is  summarized  in  Fig.  18.  Comparing  the  op¬ 
timum  configuration  at  ^»  =  6  (Fig.  14)  with  the 
optimum  configuration  at  Mach  25  (Fig,  13),  note 
that  the  Mach  25  shape  has  more  wing  sweep,  and 
pertains  to  a  con, cal  flowfield  with  a  smaller  wave 
angle,  both  of  which  are  intuitively  expected  at 
higher  Mach  number.  However,  note  from  the  flight 
conditions  listed  in  Figs.  11  and  15  that  the  body 
slenderness  ratio  at  =  6  is  constrained  to  be 
b/t  =  0.06  (analogous  to  a  supersonic  transport 
such  as  the  Concorde)  but  that  b/Z  *  0.09  is  the 
constraint  chosen  at  M»  -  25  (analogous  to  a  hydro¬ 
gen  fueled  hypersonic  aeroplane  such  as  the  British 
H0TDL).  The  two  different  slenderness  ratios  are 
chosen  on  the  basis  of  reality  for  two  different 
aircraft  with  two  different  missions  at  either  ex¬ 
treme  of  the  hypersonic  flight  spectrum.  Also  note 
in  Figs.  16-13  the  optimization  program  has  sculp¬ 
tured  a  best  optimized  configuration  with  a  spline 
down  the  center  of  tne  upper  surface--an  interest¬ 
ing  and  curious  result,  due  principally  to  the  com¬ 
peting  effects  of  minimizing  pressure  and  skin 
friction  drag,  while  meeting  the  slenderness  ratio 
constraint . 

Return  to  Fig.  15,  and  note  the  solid  symbols. 
These  pertain  to  the  values  of  C|_  and  L/D  obtained 
by  setting  the  ratio  of  specific  heats  y  to  1.1  in 
order  to  assess  possible  effects  of  high  tem¬ 
perature  chemically  reacting  flow.  The  solid  sym¬ 
bols  pertain  to  an  optimized  waverider  at  9S  =  9“ 
with  y  *  1.1.  This  is  not  necessarily  the  best 
optimum  at  Mach  25  with  y  =  1.1;  rather,  it  is  just 
a  point  calculation  to  indicate  that  high  tem¬ 
perature  effects  will  most  likely  have  a  signifi¬ 
cant  impact  on  optimized  waverider  generation,  and 
that  such  effects  are  worthy  of  future  investiga¬ 
tion.  The  detailed  aspects  of  high  temperature 
effects  are  beyond  the  scope  of  the  present  paper; 
additional  discussion  is  given  in  Ref.  7. 

As  a  final  note  in  this  section,  return  to  Fig. 
3,  and  note  that  the  solid  symbols  pertain  to  the 
present  investigation.  The  flagged  square  has  been 
discussed  earlier  as  the  data  point  for  a  half-cone 
with  a  delta  wing  at  zero  angle-of-attack ;  it  is 
not  part  of  the  present  waverider  family.  The  un¬ 
flagged  solid  square  at  M  =  4  pertains  to  a  rela¬ 
tively  large  slenderness  "ratio  of  0.087,  used  to 
generate  a  waverider  for  wind  tunnel  testing.  The 
remaining  solid  symbols,  the  circles  and  triangles, 
pertain  to  the  present  discussion.  Recall  that  the 
circles  are  for  b/t  *  0.06  (a  Concorde-like  slen¬ 
derness  ratio  for  a  low  Mach  number  configuration) , 
and  that  the  triangles  are  for  b/t  =0.09  (a  H0T0L- 
like  slenderness  ratio  for  a  high  Mach  number 
configuration).  In  the  present  section,  we  have 
discussed  results  obtained  at  Mach  6  and  25;  Fig.  3 
shows  these  plus  others  at  intermediate  Mach  num¬ 
bers.  All  of  these  cases  are  discussed  in  detail 
in  Ref.  7.  However,  in  regard  to  Fig. 3,  emphasis  is 
now  made  that  the  present  viscous  optimized  waveri¬ 
ders  produce  values  of  (L/D)  which  exceed  the  "L/D 
barrier"  discussed  in  Section  I,  and  shown  as  the 
solid  curve  in  Fig. 3.  Indeed,  the  present  waverider 

L/D  variation  is  more  closely  given  by 


shown  as  the  dashed  curve  in  Fig.  3.  Note  that 
the  two  points  given  for  M.  *  20  and  25  deviate 
away  from  the  dashed  curve.  This  is  a  Reynolds 
number  effect.  Recall  that  all  the  Mach  number- 
altitude  points  for  the  present  waveriders  are 
chosen  to  follow  a  typical  lifting  vehicle  flight 
path  through  the  atmosphere.  The  point  at  Mach 
25  is  at  very  high  altitude  (250,000  ft.),  wit£  a 
corresponding  low  Reynolds  number  (Re  -  1.4x10); 
the  flow  is  completely  laminar.  At  Mach  20,  the 
Reynolds  number  is  12  times  higher,  but  based  on 
the  transition  criterion  discussed  in  Section  II 
the  flow  is  still  completely  laminar.  Hence,  the 
laminar  skin-friction  coefficient  at  the  Mach  20 
point  in  Fig.  3  is  much  lower  (Cf  «  1//Ee)  than 
at  the  Mach  25  point,  with  an  attendant  larger 
( I/O )  at  Mach  20.  In  contrast,  the  point  at 
Mach  15  is  transitional ,  with  regions  of  both 
laminar  and  turbulent  flow,  and  hence  with 
larger  skin-friction  and  a  lower  (L/D).  In 
any  event,  the  results  given  in  Fig.  3  indicate 
that  the  present  viscous  optimized  waveriders 
produce  high  values  of  (I/O),  and  therefore  are 
worthy  of  additional  consideration  for  hypersonic 
vehicle  application. 

R.  Sensitivity  to  Transition 

Because  the  major  thrust  of  the  present  work 
is  the  inclusion  of  detailed  viscous  effects  in 
the  waverider  optimization,  the  question 
naturally  arises:  How  sensitive  are  the  present 
waveriders  to  uncertainties  in  the  location  of 
transition  from  laminar  to  turbulent  flow?  To 
address  this  question,  a  numerical  experiment  is 
carried  out  wherein  the  transition  location  was 
varied  over  a  wide  latitude,  ranging  from  all 
laminar  flow  on  one  hand,  to  almost  all  turbulent 
flow  on  the  other  hand,  with  various  cases  inbe- 
tween.  Specific  results  at  Mach  10  are  given  in 
Fig.  19;  here  values  of  (L/D)  are  given  for  opti¬ 
mized  waveriders  as  a  function  of  assumed  tran¬ 
sition  location.  The  point  corresponding  to  the 
transition  correlation  described  in  Section  II, 
0.3,  is  denoted  by  "x"  in  Fig.  19.  Other  points 
in  Fig.  19  labeled  5x,  lOx  and  15x  correspond  to 
transition  locations  that  are  5,  10  and  15  times 
the  value  predicted  by  the  transition  correla¬ 
tion.  All  the  data  given  in  Fig.  19  pertain  to 
optimized  waveriders  for  e*  *  9°,  which  yields 
the  best  optimum  at  Mach  10  for  the  usual  tran¬ 
sition  correlation.  (Note,  however,  that  8S  *  9° 
may  not  yield  the  best  optimum  for  other  tran¬ 
sition  locations;  this  effect  is  not  Investigated 
here.)  The  results  in  Fig.  19  demonstrate  a 
major  Increase  in  (l/D)  in  going  from  almost  all 
turbulent  flow  to  all  laminar  flow.  However,  for 
the  case  where  transition  is  changed  by  a  factor 
of  five,  only  a  21  change  in  L/D  results.  Even 
for  the  case  where  transition  is  changed  by  a 
factor  of  ten,  a  relatively  small  change  in  L/D 
of  lit  results.  On  the  other  hand,  the  shapes  of 
the  resulting  optimized  waveriders  are  fairly 
sensitive  to  the  transition  location,  as 
Illustrated  in  Figs.  20  and  21.  The  conclusion 
to  be  made  here  is  that  waverider  optimization 
is  Indeed  relatively  sensitive  to  transition 
location,  and  this  underscores  the  need  for 
reliable  predictions  of  transition  at  hypersonic 
speeds. 


C.  On  the  Use  of  Average  Skin  Friction  Coefficients 

The  present  detailed  viscous  analysis  computes 
the  surface  shear  stress  distributions,  and 
integrates  over  the  surface  to  obtain  the  total 
skin  friction  drag.  This  requires  a  substantial 
amount  of  computer  calculations,  and  leads  to  the 
question:  Can  an  overall  average  skin  friction 
coefficient  be  used  within  the  optimization  pro¬ 
cess  rather  than  dealing  with  the  detailed  shear 
stess  distributions?  To  address  this  question, 
consider  the  best  optimum  Mach  six  case  given  in 
Fig.  11,  which  was  orginally  calculated  with  the 
detailed  shear  stress  distributions.  From  this 
result,  an  average  skin  friction  drag  coefficient 
was  calculated  for  the  complete  configuration. 

Then  the  optimization  code  was  run  again  for  the 
same  Mach  six  case,  now  using  this  average  skin 
friction  drag  coefficient.  The  results  are  given  in 
Fig.  11  as  the  solid  symbols.  Only  a  small  dif¬ 
ference  exists  between  the  two  cases;  indeed,  the 
resulting  waverider  shapes  are  virtually  the  same, 
as  given  in  Ref.  7.  This  implies  that  if  an 
accurate  average  skin  friction  drag  coefficient  can 
be  obtained,  the  resulting  optimized  waveriders 
would  be  reasonably  valid.  However,  the  problem 
with  this  method  is  that  the  information  needed  to 
obtain  the  average  skin  friction  drag  coefficient 
is  not  known  apriori.  Moreover,  if  other  indepen¬ 
dent  means  are  used  to  obtain  an  approximate 
average  skin  friction  drag  coefficient  and  this 
approximate  average  value  is  used  in  the  optimiza¬ 
tion  process,  the  results  can  be  quite  different 
from  those  obtained  from  the  use  of  detailed  shear 
stress  distributions;  see  Ref.  7  for  more 
discussion  on  this  aspect.  This  situation,  in  com¬ 
bination  with  the  sensitivity  to  transition 
demonstrated  in  the  previous  section,  seems  to  dic¬ 
tate  the  necessity  of  using  the  detailed  shear 
stress  distributions  rather  than  some  approximate 
average  value  of  skin  friction  drag  coefficient  for 
obtaining  the  proper  optimized  waveriders. 


As  a  final  note,  it  is  interesting  to  pose  the 
question:  if  the  skin  friction  is  deleted  from  the 
present  analysis,  what  type  of  optimized  inviscid 
waverider  configuration,  with  a  constraint  on  sTen- 
derness  ratio.  Is  produced?  To  examine  this 
question,  the  present  computer  code  was  run  without 
skin  friction  as  part  of  the  optimization  process, 
covering  the  range  of  Mach  number  from  6  to  25.  A 
typical  result  for  the  inviscid  optimized  configu¬ 
ration  is  shown  in  Fig. 22.  Here  we  see  essentially 
a  wedge-like  caret  waverider,  such  as  the  classic 
configuration  generated  by  the  two-dimensional  flow 
behind  a  planar  oblique  shock  wave,  as  discussed  by 
Nonweiler  ’  .  This  clearly  Indicates  that  the 
optimized  inviscid  waverider  with  slenderness  ratio 
as  the  constraint  is  Indeed  a  caret  wing.  The 
result  shown  in  Fig.  22  is  produced  by  the  present 
conical  flow  analysis  as  a  "limiting  case",  wherein 
the  optimum  shape  Is  seeking  the  flattest  portion 
of  the  conical  shock  wave.  To  see  this  more 
clearly,  return  to  Fig.  6.  The  resulting  inviscid 
waveriders  are  being  generated  by  relatively  flat 
streamsurfaces  at  the  extreme  back  and  bottom  of 
the  generating  conical  flow-field  —  where  the 
shock  radius  of  curvature  Is  the  largest  and  the 
flow  is  closest  to  being  two-dimensional. 
Consequently,  the  inviscid  configurations  are  tiny 
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shapes  compared  to  the  scale 
of  the  flowfleld  in  Fig.  6,  and  they  are 
“squeezed"  Into  a  tiny  area  at  the  bottom  of  the 
shock  base.  In  turn,  due  to  the  logic  of  the 
existing  conical  flow  code,  only  a  few  pressure 
and  shear  stress  points  are  calculated  on  the 
surface  of  these  tiny  waveriders,  raising 
questions  about  the  numerical  accuracy  of  the 
calculation  of  their  lift  and  drag.  Therefore, 
no  further  discussion  about  the  inviscid  opti¬ 
mized  waveriders  will  be  given  here,  except  to 
emphasize  again  that  a  two-dimensional  caret  wing 
seems  to  be  the  optimum  1 nvi scid  shape  that  is 
predicted  by  the  present  conical  flow  analysis. 

IV  CONCLUSIONS 

In  comparison  to  previous  optimized  waverider 
analyses,  the  present  work  is  the  first  to 
Include  detailed  viscous  effects  within  the  opti¬ 
mization  process.  From  this  work,  the  following 
major  conclusions  are  made: 

1.  The  resulting  family  of  viscous  hyper¬ 
sonic  waveriders  yields  predicted  high  values 
of  (L/o )  which  break  the  "L/0  barrier" 
discussed  in  Section  I. 

2.  The  optimization  process  for  the  viscous 
waveriders  results  in  distinctly  different 
shapes  compared  to  previous  work  with 
inviscid-designed  waveriders. 

3.  The  fine  details  of  the  viscous  solution, 
such  as  how  the  shear  stress  is  distributed 
over  the  surface,  and  the  location  of  tran¬ 
sition,  are  crucial  to  the  details  of  the 
resulting  waverider  geometry. 
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FIG.  2:  Performance  comparison  of  various 
generic  hypersonic  configurations 
(  from  Ref.  3). 
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FIG  3:  Maximum  llft-to-drag  ratio  comparison 
for  various  hypersonic  configurations. 


FIG  1:  Newtonian  results  for  a  flat  plate. 


IMCRItIC 


\\\ 

V'V 


Pionor.aiteched 

shock  wove 


Schematic  of  a  simple  caret  wing 
waverider. 
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FIG.  5:  Spherical  coordinate  system  for  cone 
flow  calculations,  and  cartesian 
coordinate  system  for  waverider 
design. 
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FIG  6:  Perspective  view  of  a  conical 
waverider  and  generating  shock. 
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FIG.  7:  Front  view  of  a  conical  waverider 
with  the  generating  cone  and  shock. 


FIG  8:  Expansion  cylinder  and  expansion 
domain  bounded  by  the  Mach  cone. 
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FIG.  9:  Example  of  the  optimization  process 
for  a  function  of  two  variables. 
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FIG.  15:  Results  For  a  series  of  optimized 
waver iders  at  Mach  25. 
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FIG.  16:  Front  views  of  a  series  of  optimized  waveriders  at  Mach  25. 
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FIG.  19: 


Lift-to-drag  comparison  of  optimized 
Mach  10  waveriders  designed  with  various 
boundary  layer  transition  criteria. 
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Front  views  of  optimized  Mach  10  waveriders  designed  with  various 
boundary  layer  transition  criteria. 
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III.  SHOCK  WAVE/VORTEX  INTERACTION 


This  research  is  the  Ph.O.  dissertation  topic  for  Griffin  Corpening, 
scheduled  for  completion  in  May,  1988.  The  problem  involves  the  study  of  the 
interaction  between  an  oblique  shock  wave  and  a  vortex  at  hypersonic  speeds. 

This  is  modeled  as  a  three-dimensional  inviscid  flow.  The  solution  technique 
involves  the  appliecation  of  computational  fluid  dynmaics;  specifically,  a 
time-marching,  finite-volume  approach  is  used,  including  upwind  differencing. 

The  computati onal  method  is  patterned  after  that  of  Peter  Gnoffo  at  the  NASA 
Langley  Research  Center,  and  has  been  throughly  tested  for  two-dimensional 
normal  and  oblique  shocks,  intersecting  shocks,  and  slip  lines  intersecting 
oblique  shocks.  The  three-dimensional  method  has  not  been  tested  as  throughly 
due  to  computer  time  and  memory  constraints,  but  tests  performed  so  far  indicate 
the  model  to  be  operating  properly. 

In  addition  to  the  actual  code  development  itself,  it  was  necessary  to  write 
a  number  of  programs  to  enable  the  flow  field  to  be  examined  graphically.  Due 
to  the  three-dimensional ity ,  this  was  no  trivial  undertaking  (on  the  order  of 
three  to  four  months  was  required  to  do  this).  The  flow  field  can  now  be 
displayed  by  two-dimensional  slices  perpendicular  to  any  of  the  three  axes  and 
at  various  locational  stations  along  the  axis.  The  two-dimensional  slice  can  be 
reduced  to  contour  plots,  surface  plots,  velocity  vector  or  streamline  plots. 
Curves  can  also  be  plotted  showing  the  velocity,  density,  etc.  change  through 
the  slice.  In  addition,  three-dimension  streamlines  can  be  calculated  and 
plotted.  All  this  combines  to  give  a  very  flexible  and  complete  way  to 
graphically  display  the  flowfield. 

Some  sample  results  are  given  as  follows.  (The  figures  are  in  a  preliminary 
form  which  will  be  formalized  for  the  dissertation. )  The  following  figures 
depict  a  Mach  5  flowfield  in  which  a  vortex  is  inte-acting  with  a  shock  wave. 

The  flowfield  is  depicted  by  showing  a  series  of  contour  plots  of  density. 
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These  plots  are  taken  along  planes  which  are  perpendicular  to  the  3  axes.  A 
Three-Dimensional  Perspective  figure  is  also  included  to  help  the  reader  orient 
the  contour  plots  to  the  flowfield. 


The  freestream  flow  direction  is  in  the  positive  X  direction.  The  shock 
wave  starts  along  the  Z-axis  at  Y  equal  to  0  and  proceeds  up  into  the  flow  at  an 
angle  of  46.4  degrees  to  the  X-Z  plane  and  proceeds  downstream. 
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IV.  CONCLUSION 

The  Fellowship  Program  in  Hypersonic  Aerodynamics  sponsored  at  the 
University  of  Maryland  by  ARO  has  been  very  productive,  far  beyond  its  initial 
scope.  The  work  of  two  exceptional  Ph.D.  students  is  summarized  in  this  report 
A  similar  fellowship  program  (now  involving  one  Ph.D.  student)  is  currently  in 
progress,  and  will  be  reported  at  the  proper  time. 


